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I.  INTRODUCTION 

Continuing  safe  life  service  of  infrastructure  is  a  major  challenge  in  current  society  and 
culture.  An  immediate  focus  on  aging  aircraft  has  generated  a  multitude  of  problem  analyses  and 
resolution  approaches.  The  development  of  fracture  mechanics  as  an  analysis  tool  has  enabled 
quantification  of  scheduled  aircraft  inspection  and  maintenance  and  the  quantification  damage 
tolerance  characteristics  and  service  margins  for  crack-like  anomalies.  Crack-like  anomalies  due 
to  variance  in  material  quality,  service  usage  (fatigue  and  physical  damage)  and  environmental 
initiation  are  the  primary  initial  consideration  in  assessing  and  quantifying  structural  integrity  of 
engineering  materials,  structures,  components  and  systems.  Crack-like  anomalies  are  the  cause 
for  approximately  80%  of  the  aircraft  that  are  retired  from  service  due  to  degradation  in  structural 
integrity.  Major  technology  initiatives  have  focused  on  detection  and  quantification  of  crack-line 
anomalies  with  resultant  high  confidence  in  aircraft  life  extension  and  major  cost  savings  in 
operational  readiness  and  aircraft  replacement. 

Corrosion  is  frequently  cited  as  a  major  factor  in  aircraft  maintenance,  component 
replacement,  and  life-cycle  cost.  Indeed  corrosion  is  a  source  of  crack  initiation  with  resultant 
degradation  in  structural  integrity  and  is  of  concern  due  to  loss  of  material  and  hence  general 
degradation  in  structural  integrity.  The  engineering  parameter  most  often  used  for  assessing 
degradation  due  to  corrosion  is  loss  of  material  /  thinning.  Fortunately,  load  redistribution  and 
general  margins  provide  damage  tolerance  due  to  corrosion  over  a  wide  range  of  corrosion 
conditions.  Unfortunately,  corrosion  damage  is  often  hidden  and  is  not  detected  in  a  timely 
manner  to  enable  optimum  scheduled  arrest,  maintenance  and  replacement.  Continuing  life 
extension  of  aging  aircraft  has  increased  the  need  for  both  timely  detection,  quantification  and 
analysis  of  corrosion  (visible  and  hidden)  for  safe-life  aircraft  systems  management. 

Development  and  evolution  of  a  measurement  technology  is  often  based  on  what  can  be 
measured  and  on  criteria  described  in  terms  of  what  can  be  measured.  In  the  case  of  corrosion, 
wall  thinning  in  an  aircraft  skin  can  be  measured  using  available  ultrasonic,  eddy  current  and 
radiometric  tools.  A  loss  in  wall  thickness  can  then  be  related  to  an  increase  in  a  local  stress 
concentration  using  load  /  stress  per  unit  area  analyses.  Unfortunately,  corrosion  also  occurs  at 
interfaces  and  in  the  second  layer  of  an  aircraft  structure  that  may  not  be  accessible  and  is  not 
reliability  assessed  by  conventional  nondestructive  evaluation  tools.  The  problem  is  further 
complicated  by  the  lack  of  structural  engineering  analysis  tools  and  acceptance  criteria.  Reliable 
nondestructive  measurement  /  assessment  tools  are  necessary  for  economical  management  of 
corrosion  damage  and  quantified  engineering  acceptance  criteria  are  necessary  to  realize  benefit 
from  early  detection,  measurement  and  characterization.  This  “chicken  and  egg”  dilemma  has 
been  approached  by  a  combination  of  exploratory  approaches  to  corrosion  detection  and 
measurement  and  by  exploratory  approaches  to  understanding  and  quantifying  the  effects  of 
corrosion  damage  on  structural  integrity.  The  work  described  herein  addresses  elements  of  an 
exploratory  approach  to  reliable  nondestructive  measurement  of  corrosion  damage  characteristics 
and  a  cursory  assessment  of  related  work  in  establishing  acceptance  criteria. 
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The  nondestructive  measurement  work  is  an  extension  of  work  initiated  by  the  late  John 
Moulder  [REF  1]  in  the  application  of  pulsed  eddy  current  measurements  to  a  two-layer  structure 
that  has  been  modified  by  corrosion  damage.  The  Moulder  work  was  the  result  of  experimental 
observations  and  a  deep  knowledge  of  practical  eddy  current  interactions  and  applications.  The 
advancements  provided  by  the  current  project  were  to  develop  and  measurement  model  for 
predicting  detection  capability  and  improvements  in  the  “probability  of  detection”  analysis  tools 
used  in  assessing  the  capability  and  reliability  of  nondestructive  measurements  and  applications. 
Finally,  an  assessment  of  potential  for  commercial  implementation  of  the  measurement  technique 
and  instrumentation  was  initiated. 

The  project  reported  herein  was  divided  into  four  tasks: 

•  Assessment  and  definition  of  the  need  for  corrosion  measurement  based  on  the  pulsed 
eddy  current  method; 

•  Development  and  validation  of  a  measurement  model  using  the  pulsed  eddy  current 
method; 

•  Improvements  in  the  analysis  code  used  for  the  assessment  of  nondestructive 
measurement  capability  and  reliability  (probability  of  detection);  and 

•  Assessment  and  planning  commercial  implementation  of  the  pulsed  eddy  current  method 
(Planned  Phase  II) 

II  SIGNIFICANCE  OF  THE  HIDDEN  CORROSION  PROBLEM  AND 

OPPORTUNITY  FOR  IMPLEMENTATION  OF  A  PULSED  EDDY  CURRENT 
MEASUREMENT 

A.  GENERAL 

The  NRC  Report,  “Aging  of  US  Air  Force  Aircraft”  [REF  2]  identifies  the  most 
important  NDE  needs  as  (1)  “Damage  detection  of  fatigue  cracks  under  fasteners”;  (2)  “Hidden 
Corrosion”;  (3)  “Detection  of  small  cracks  associated  with  wide  spread  fatigue  damage  (WFD); 
and  (4)  “Corrosion  in  multi-layer  structures”.  Although  the  need  has  been  broadly  identified, 
specific  identification  and  characterization  of  that  need  are  necessary  to  assess  the  primary 
factors  involved  in  establishing  requirements  for  measurement  and  matching  those  requirements 
to  available  technology.  Commercialization  of  a  measurement  capability  (tool,  system  or  service) 
must  include  analysis  of  the  need,  practicality  and  economics  of  a  broad  application. 

Existing  ultrasonic,  eddy  current  and  radiometric  methods  adequately  measure  wall 
thinning  in  the  first  layer  of  an  aircraft  structure.  In  addition,  wide  area  imaging  by  pulsed 
infrared  methods  provide  a  semi  quantitative  measure  of  corrosion  patterns  and  location. 
Unfortunately,  there  is  no  standardized  acceptance  criteria  or  methods  of  using  the  data  provided 
by  these  measurements.  Indeed,  each  aircraft  manufacturer  and  user  has  an  individualized 
method  of  assessing  and  using  the  information  provided.  A  part  of  the  lack  of  standardization  is 
due  to  the  variation  in  use  and  need  between  commercial,  private  and  military  aircraft.  The  same 
difficulty  is  evident  in  identifying  criteria  for  detection  and  quantification  of  hidden  corrosion  in 
all  aircraft  structures.  The  structures  engineering  analysis  task  is  difficult  and  methods  of  analysis 
are  immature.  Indeed,  current  plans  for  structures  assessment  (with  respect  to  hidden  corrosion) 
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include  a  detailed  tear  down  and  analysis  of  retired  aircraft  to  try  to  understand  the  problem  and 
to  develop  criteria  and  analytical  methods.  Two  parameters  that  are  evident  in  analysis  needs  are 
the  gap  between  interfaces  in  multi-layer  aircraft  skin  and  corrosion  damage  in  the  second  layer. 
The  pulsed  eddy  current  measurement  method  is  currently  the  only  method  identified  which  may 
meet  the  needs  for  identification  and  quantification  of  characteristic  damage  conditions. 

B.  NEED  AND  ACCEPTANCE  CRITERIA 

Survey  of  industry  aircraft  structure  engineering  experts  reaffirmed  the  lack  of  quantified 
acceptance  criteria  for  second  layer  corrosion  and  developmental  state  of  understanding  and  tools 
for  analysis  and  management  of  corrosion  damage.  The  most  promising  program  involves 
disassembly  of  a  retired  commercial  aircraft;  visual  analysis  and  photographic  recording  of 
conditions  found;  sub-element  testing  of  selected  components  to  assess  residual  strength  and 
fatigue  life;  and  an  exploratory  analysis  to  relate  corrosion  pitting  patterns  to  residual  fracture 
toughness.  Near  term  closure  is  not  expected. 

In  summary,  there  is  a  recognized  need  to  detect  and  quantify  progressive  corrosion 
damage  in  aircraft  structures.  Current  practices  rely  on  deterministic  engineering  judgment  for 
repair  /  replacement  of  damaged  areas.  In  addition,  corrosion  damage  to  military  aircraft  differs 
significantly  from  that  experienced  in  commercial  and  private  aircraft.  The  differences  are 
primarily  due  to  differences  in  materials  and  alloys,  loading  /  flight  profiles,  design  margins,  base 
facilities  environments  and  maintenance  practices. 

C.  DEFACTO  CRITERIA  BASED  ON  INSPECTION  /  DETECTION  CAPABILITIES 

In  the  absence  of  mature  analytical  methods  and  quantified  acceptance  criteria,  the 
capability  of  available  inspection  methods  has  traditionally  been  the  defacto  acceptance  criteria. 
Use  of  such  tools  aids  in  management  of  engineering  structures  and  maintenance  operations  and 
supports  increasing  understanding  of  the  conditions  and  effects  of  the  anomaly  being  assessed. 
Pulsed  eddy  current  measurement  methods  have  been  demonstrated  to  be  capable  of: 

•  Detecting  corrosion  damage  in  the  form  of  thinning  in  the  first  layer  of  an  aircraft 
structure; 

•  Detecting  gaps  between  the  first  and  second  layer  of  an  aircraft  structure  (note 
that  the  gap  is  often  increased  by  the  accumulation  of  corrosion  products  at  the 
interface);  and 

•  Detecting  corrosion  damage  in  the  form  of  thinning  in  the  second  layer  of  an 
aircraft  structure. 

Aircraft  owners  /  operators  will  be  reluctant  to  invest  in  instrumentation  to  detect  and 
measure  second  layer  corrosion  until  a  generally  accepted  acceptance  criteria  is  developed  and 
economic  methods  of  corrosion  arrest  and  remedial  action  are  developed.  Technology  transfer  of 
pulsed  eddy  current  measurement  methods  to  industry  will  therefore  depend  on  its  capabilities 
for  providing  useful  measurements  in  detecting  extreme  conditions  and  its  usefulness  in 
understanding  the  conditions  and  effects  of  corrosion  damage  on  structural  integrity.  The  primary 
objectives  of  the  program  reported  herein  were  therefore  to: 
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•  Increase  understanding  of  the  capabilities  of  the  pulsed  eddy  current  method  in 
materials  evaluation  as  applied  to  corrosion  damage; 

•  Increase  understanding  of  the  pulse  eddy  current  method  with  the  objective  of 
establishing  confidence  in  its  use  as  a  measurement  method; 

•  Increase  awareness  of  the  method  as  a  useful  tool  for  nondestructive  measurement  and 
assessment;  and 

•  To  explore  technology  transfer  potential  including  the  development  and  marketing  of 
a  commercial  instrument  for  general  use  in  the  public  domain. 

III.  THE  PULSED  EDDY  CURRENT  (PEC)  METHOD 

Although  several  NDE  methods  are  available  to  detect  and  assess  corrosion  damage  in 
the  first  layer  of  aircraft  structures,  the  pulsed  eddy  current  (PEC)  method  developed  by  Center 
for  Nondestructive  Evaluation  (CNDE)  at  Iowa  State  University  (ISU)  is  the  only  method  that 
has  been  demonstrated  to  be  capable  of  quantifying  interface  gap  (varying  amounts  of  corrosion 
products)  or  reliable  detection  in  the  second  layer  [REF  1].  Second  layer  and  hidden  corrosion 
have  been  identified  as  major  aging  aircraft  problem  conditions.  Proof  of  principle  and 
applicability  of  the  PEC  method  been  demonstrated  and  validated  in  related  work. 

Pulsed  eddy  current  work  at  CNDE  was  initiated  by  the  late  John  Moulder  and  resulted  in 
the  development  of  a  novel  approach  to  application  and  analysis  of  the  technology.  The  ISU 
work  is  highly  leveraged  by  combined  support  by  this  project,  by  the  Federal  Aviation 
Administration  (FAA)  and  by  private  industry.  Joint  benefits  have  been  realized  by  concurrent 
work  on  different  elements  of  pulsed  eddy  current  research,  methods  validation  and 
demonstration  on  industry  products. 

Pulsed  eddy  current  measurement  is  not  a  new  technology.  Modem  integrated  circuit 
components  and  digital  design  capabilities  have  greatly  increased  measurement  instrument 
stability  and  options  for  complex  measurements.  The  parameters  measured  and  the  causal  model 
relationships  developed  in  our  work  are  new  innovations  that  have  enabled  assessment  of 
corrosion  conditions.  To  our  knowledge,  the  method  is  a  unique  advancement  in  eddy  current 
technology  and  experimental  assessments  have  shown  detection  capabilities  for  hidden  corrosion 
and  corrosion  a  second  layers  in  the  presence  of  a  varying  interface  gap.  The  PEC  technology  is 
an  advanced  form  of  the  usual  continuous-wave  (CW)  eddy  current  (EC)  method,  which  is  the 
traditional  choice  of  NDE  technique  for  many  inspection  applied  to  aircraft  structures.  However, 
in  order  for  CW  EC  instruments  to  qualify  for  the  above  types  of  critical  inspections,  they  must 
operate  at  low  frequencies  for  deep  penetration  and  at  multiple  frequencies  for  depth  resolution. 
These  operational  requirements  are  both  time-consuming  and  inconvenient.  The  time-domain  EC 
method  provides  an  excellent  resolution  of  these  shortcomings  because,  by  nature,  it  can  acquire 
the  equivalent  of  low-frequency  and  multiple-frequency  data  in  a  single  pulse.  Being  a  time- 
domain  method,  pulsed  eddy  current  provides  opportunities  to  employ  time-gating,  which  has 
shown  excellent  discrimination  capabilities  against  interfering  signals  from  fasteners  and 
subsurface  structures. 
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The  state  of  the  art  of  the  ISU  PLC  technology  is  found  in  [REF  1,6].  This  work 
includes  the  Cheng-Dodd-Deeds  (CDD)  theory  behind  the  model-based  corrosion 
characterization.  Pulsed  eddy-current  systems  have  important  advantages  from  the  operator 
perspective.  They  can  be  constructed  from  simple,  relatively  inexpensive  equipment.  Since  the 
method  is  based  on  a  familiar  technology,  problems  with  implementation  of  the  technique  in  a 
USAF  overhaul  or  field  level  NDE  facility  by  USAF  personnel  would  be  minimized.  By  adding  a 
scanning  capability  to  the  instrument,  the  results  can  be  displayed  in  an  image  format,  which  is 
much  easier  to  interpret  and  understand.  Since  pulsed  signals  contain  inherently  a  broad  range  of 
frequencies,  pulsed  eddy  current  methods  are  akin  to  swept-frequency  measurement  methods. 

Yet  the  pulsed  eddy-current  technique  is  much  faster:  100  milliseconds  to  acquire  a  pulsed  eddy- 
current  signal  compared  to  several  minutes  to  complete  a  swept-frequency  measurement. 
Operation  of  a  prototype  instrument  has  been  demonstrated  at  Tinker,  Kelly,  McClellan,  and 
Robins  Air  Force  Bases.  It  was  singled  out  as  one  of  the  most  promising  new  technologies  for 
corrosion  detection  in  the  second  round  of  blind  trials  conducted  by  ARINC,  Inc.  for  Tinker  AFB 
[REF  3,4,5], 

The  original  PEC  method  demonstrated  used  a  bobbin  coil  for  both  generation  and 
detection.  A  false  color  image  of  the  test  object  was  generated  from  the  time-gated  output  signal 
produced  as  shown  in  Figure  1 .  Figures  2  and  3  describe  various  features  of  system  functionality. 
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Figure  1.  Pulsed  eddy  current  instrument  control  and  display  panels  (software).  The  difference  between  the 
return  and  reference  signals  is  plotted  as  a  function  of  time  in  the  A-scan  display.  For  2D  scans,  the  instrument 
displays  a  C-scan  image  based  on  the  peak  values.  The  software  allows  users  to  display  the  time-gated  C-scan  image 
for  a  selected  range  of  zero-crossing  times. 
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Figure  2.  Plots  to  show  the  usage  of  the  time-of-flight  information,  which  can  determine  the 
severity  of  corrosion  as  well  as  the  relative  location  in  the  multi  layer. 
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Figure  3.  Time-gated  vs.  non-gated  C-scan  images.  These  plots  demonstrate  that  the  time-gating 
method  effectively  discriminates  corrosion  signals  from  other  signals  of  geometrical 
origins,  such  as  rivets  and  holes.  [REF  1], 
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Instrument  improvements  were  made  during  the  current  program  on  a  companion  project. 
Detector  coils  were  replaces  by  Hall-  sensors  to  improve  depth  of  penetration  and  remove 
sensitivity  to  sensor(primarily  temperature)  changes.  The  detector  board  was  upgraded  to 
accommodate  the  Hall-sensor  and  to  improve  signal  outputs.  Detector  software  was  upgraded  to 
Win-32  and  analysis  /  processing  improvements  were  made  to  remove  stringer  and  tear-strap 
signals  and  to  improve  data  storage.  Resultant  system  performance  was  improved  as  shown  in 
Figure  4.  The  depth  of  penetration  was  increased  using  the  Hall-sensor  as  shown  in  Figure  5. 


•  Detection  capability:  less  than  4%  material  loss  in  Vi"  Al  plate. 

*  Depth  penetration  can  be  increased,  at  the  expense  of  measurement 
time,  using  a  larger  probe. 


Figure  4.  Improved  detection  capability  using  the  Hall-Sensors. 


Slot  length:  0.5” 

Slot  width:  0.015” 

Slot  plate  thickness:  0.025” 
Slot  plate 

conductivity:  34.5%  IACS 


Figure  5  Improvements  in  depth  of  penetration  using  the  Hall-sensor 
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A  summary  of  the  concurrent  hardware  and  experimental  tasks  completed  under  funding 
from  the  FAA  is  included  for  reference  in  APPENDIX  A. 


IV.  MEASUREMENT  MODEL  DEVELOPMENT  FOR  PULSED  EDDY  CURRENT 
APPLICATIONS 

A  major  contribution  and  output  of  the  current  project  was  the  development  of  a 
measurement  model  for  pulsed  eddy  current  application.  The  model  is  a  useful  tool  in 
understanding  various  aspects  of  pulsed  eddy  current  technology  and  in  quantifying  expected 
results  for  various  applications.  The  predictive  model  is  also  a  useful  input  to  prediction  of 
detection  capabilities  and  reliability  using  the  probability  of  detection  (POD)  models.  The  model 
predicts  output  signal  levels  in  various  applications  and  must  be  used  in  combination  with 
experimentally  determined  noise  measurements  for  use  in  the  POD  models.  Noise  levels  are 
input  for  laboratory  and  field  measurements  and  provide  the  anchor  for  predicting  detection 
capabilities  in  the  laboratory  /  field  applications  [REF  7,8,9,]. 

John  R.  Bowler  completed  the  measurement  model  as  an  independent  task  as  a  part  of 
the  CNDE  support  work.  The  stand-alone  task  report  for  this  work  is  included  in  Appendix  B. 

The  Software  program  for  this  model  is  available  from  the  Center  for  Nondestructive  Evaluation 
at  Iowa  State  University.  [REF  10].  The  output  of  the  software  model  is  data  in  the  form  shown 
in  Figure  5.  Instructions  for  use  of  the  software  model  are  described  in  APPENDIX  C. 

V.  PROBABILITY  OF  DETECTION  (POD)  MODEL  INTERFACE  AND  UPGRADE 

Probability  of  Detection  (POD)  is  the  standard  for  quantifying  the  capability  and 
reliability  of  various  nondestructive  evaluation  (NDE)  procedures  and  applications.  The  original 
work  in  development  of  the  method  was  completed  using  a  moving  average  method  of 
quantifying  the  POD  function.  This  method  requires  a  large  amount  of  data  and  is  uneconomical 
for  most  applications.  The  form  of  the  function  is  described  as  the  merger  of  the  distributions  of 
signal  and  noise  in  a  practical  application.  Signal  and  noises  distributions  are  expected  to  be 
Gausian  in  form,  thus  the  merger  of  the  two  distributions  can  be  described  in  the  form  a  log 
logistic  model.  This  method  was  proposed  by  Berens  [REF  11]  and  has  become  the  standard  for 
POD  data  analysis. 

The  POD  (detection)  model  is  based  on  an  assumption  that  repetitive  measurements  of  a 
single  flaw  and  the  associated  measurement  noise  are  normally  distributed;  and  that  the  measured 
values  provide  a  means  for  discriminating  between  background  measurement  noise  and  a  flaw 
signal  [REF  7,12,13],  Figure  6,  shows  a  typical  desired  response  with  clear  discrimination 
between  signal  and  noise.  The  measurement  distributions  are  the  result  of  repetitive 
measurements  made  on  a  flaw  of  a  single  size.  A  number  of  measurements  are  required  to 
establish  confidence  in  the  repeatability  of  the  measurement  process.  The  corresponding 
confidence  level  and  required  sample  size  are  taken  directly  from  statistical  sampling  tables  and 
are  based  on  random  sampling  (for  example,  MIL  HDBK  5). 
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As  the  flaw  size  decreases,  the  measured  signal  magnitude  decreases  and  approaches  the 
noise  measurement  levels  as  shown  in  Figure  7.  Discrimination  between  signal  and  noise  are 
more  difficult  and  “FALSE  CALLS”  will  be  made  when  the  noise  signal  is  interpreted  as  a  flaw 
signal.  If  the  acceptance  threshold  is  set  low  enough  (point  A)  to  assure  detection  of  all  flaws,  a 
significant  number  of  “FALSE  CALLS”  will  result  due  to  noise  being  interpreted  as  a  FLAW 
signal  level.  Conversely,  if  the  acceptance  threshold  is  set  high  enough  to  avoid  “FALSE 
CALLS”,  some  “FLAWS”  will  be  missed  and  the  probability  of  detection  will  be  low.  The  useful 
threshold  discrimination  point  (specificity)  for  a  measurement  method  is  that  point  at  which  the 
probability  of  detection  (POD)  exceeds  a  minimum  value  (by  convention  the  90%  point)  with  no 
false  calls  (or  an  accepted  level  of  “false  calls”).  Some  flaws  that  are  smaller  than  the  established 
threshold  POD  value  will  be  detected  due  to  flaw-to-flaw  variances,  but  the  measurement  method 
will  not  reliably  detect  flaws  that  are  below  the  POD  value.  Validation  of  the  capability  of  an 
NDE  procedure  requires  characterization  of  the  signal  and  noise  response  values  and  the 
establishment  of  a  useable  discrimination  threshold,  (acceptance  criteria  for  the  procedure).  The 
decision  threshold  level  is  most  often  established  by  comparative  measurement  to  a  reference 
artifacts  (“calibration”)  and  is  the  basis  for  reproducing  the  measurement  parameters  for  an  NDE 
procedure.  The  causal  relationship  response  of  the  measurement  to  variation  in  flaw  size  is 
established  experimentally  by  repetitive  measurements  of  flaws  of  varying  size  (ideally  the 
relationship  is  linear  or  log-linear). 


SIGNAL  MAGNITUDE 
(AMPLITUDE) 


Figure  6.  Relationship  between  measurement  signal  and  background  /  measurement  noise 
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SIGNAL  MAGNITUDE 
(AMPLITUDE) 

Figure  7.  Signal  and  noise  overlap  at  a  small  flaw  size 

The  standard  statistical  model  for  describing  two  overlapping  Gausian  distributions  is  the  log- 
logistics  (log-odds)  model  that  is: 


PoD(«)  =  exp«x+ftlnaL  or  ln 
l+exp(a+plna) 
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Some  investigators  with  similar  results  also  use  the  lognormal  or  Weibull  model.  Both 
POD  methods  assume  a  log  linear  relationship  between  increasing  signal  level  and  increasing 
flaw  size.  This  relationship  is  valid  for  most  NDE  modalities  when  the  flaw  is  large  with  respect 
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to  the  interrogating  probe  and  with  those  modalities  that  depend  on  visual  detection  such  as 
visual  and  penetrant  methods. 

An  integrated  NDE  measurement  model  (such  as  that  developed  for  the  pulsed  eddy 
current  application)  may  be  used  to  predict  output  signal  levels  under  varying  instrument  and 
operational  parameters  and  may  therefore  be  linked  to  predict  the  probability  of  detection  (POD) 
under  varying  instrument  /  operational  conditions.  The  measurement  model  must  also  account  for 
“noise”  signals  such  that  the  threshold  discrimination  level  (acceptance  criteria)  can  be  set  to 
detect  all  measured  anomalies  without  reporting  false  calls.  The  “noise”  component  may  be 
input  from  models  or  from  experimental  measurement  data  [REF  7,8,9].  An  acceptance  signal 
threshold  must  then  be  established  to  provide  the  basis  for  “detection”  or  a’  miss  in  the  case  of  a 
‘Hit”  /  “Miss”  analysis  or  as  a  boundary  condition  for  analysis  of  discrete  signal  data.  The  linked 
modeling  sequences  constitute  a  predictive  POD  modeling  tool.  An  interim  step  may  be 
necessary  to  transform  the  predicted  output  to  the  log  linear  form  that  is  necessary  for  the 
logistics  or  lognormal  POD  models.  Attention  to  the  transform  is  especially  required  for 
detection  of  small  flaws  and  for  detection  involving  special  boundary  conditions,  such  as  an  edge 
at  a  corrosion  damage  site. 

Basic  POD  software  (log-logistic  model)  was  modified  in  Microsoft  Excel™  Macro 
[REF  1 3  ]  form  to  provide  an  engineering  user  tool  for  analysis  of  “Hit  /  Miss”  data.  The  basic 
software  was  used  as  the  basis  for  analyses  of  all  data  in  the  “NDE  Capabilities  Data  Book” 

[REF  12].  Since  different  versions  of  Microsoft  Excel™  require  different  procedures  for  plotting, 
specific  instructions  for  three  different  versions  were  required.  The  Microsoft  Excel™  Macro 
method  was  selected  for  this  software  package  to  provide  direct  interface  with  most  data  in 
tabular  /  spread  sheet  form.  This  provides  a  convenient  interface  to  the  eddy  current  measurement 
model  data  and  to  experimental  data.  The  POD  Excel  Macro  is  available  from  D&W  Enterprises 
[REF  14]. 

The  output  of  the  new  Macro  analysis  software  is  shown  in  Figure  8,  where  the  threshold 
Snde  value  (by  convention)  is  that  point  at  which  the  curve  crosses  the  90-percentile  level.  The 
output  results  are  identical  to  that  used  in  the  “NDE  Capabilities  Data  Book”.  The  analysis  and 
plotting  are  considerably  faster  than  then  the  basic  model  and  are  now  able  to  interface  with 
various  version  of  Microsoft  Excel™.  Instructions  for  use  of  the  POD  Macro  are  documented  in 
APPENDIX  D.  Plots  of  data  from  the  pulsed  eddy  current  model  were  not  completed  due  to  the 
requirement  for  experimentally  measured  noise  data.  For  purposes  of  this  initial  assessment,  the 
noise  level  is  assumed  to  be  at  the  baseline  for  the  pulsed  eddy  current  model  output. 
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Figure  8.  A  typical  POD  curve  output  from  the  POD  Macro 

VI.  TECHNOLOGY  TRANSFER  ANALYSIS  AND  SUMMARY 

A  major  objective  of  the  project  was  to  assess  marketability  and  transfer  of  the 
technology  to  the  commercial  sector.  Several  NDE  instrument  producers  expressed  interest  in 
producing  an  instrument  under  contract  as  a  research  instrument  or  to  customer  specific 
requirements.  Such  production  would  include  the  cost  of  instrument  design  specific  to  Air  Force 
needs  and  provisions  for  minimum  quantity  buy  would  be  included.  This  is  typical  of  Air  Force 
specific  NDE  instrumentation  that  are  used  in  other  inspection  applications  and  are  used  in 
general  Air  Force  maintenance  operations.  The  same  instrument  manufacturers  would  consider 
licensing  and  production  of  a  commercial  instrument  at  such  time  that  firm  requirements  are 
generated  by  the  aircraft  (original  equipment  manufacturers  OEM’s)  industry  and  of  profitable 
market  can  be  projected.  Although  the  technology  uniquely  addresses  an  Air  Force  need, 
commercialization  is  not  recommended  at  this  time.  A  Phase  II  program  will  therefore  be  delayed 
to  allow  for  further  development  of  the  technology  and  equally  important,  to  allow  for 
development  of  engineering  acceptance  criteria  and  analysis  methods  to  utilize  the  measurements 
made.  Requirements  for  measurement  will  not  be  forthcoming  until  engineering  acceptance 
criteria  and  analysis  methods  are  developed  and  generally  accepted  by  Air  Force  Maintenance,  by 
aircraft  manufacturers,  and  by  airlines.  The  pulsed  eddy  current  measurement  technology  is  still 
evolving  as  evidenced  by  the  replacement  of  the  coil  with  a  Hall  sensor  detector,  by  the  industry 
evaluation  /  assessment  program  with  Boeing  [REF  16],  and  by  the  maturity  of  the  pulsed  eddy 
current  instrument  design  and  acceptance  test.  Major  progress  in  maturation  of  the  technology 
was  provided  in  the  form  of  the  measurement  model  and  the  POD  macro  programs  that  are 
necessary  for  analysis. 
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In  summary,  the  pulsed  eddy  current  method  offers  a  unique  capability  for  detection  and 
quantification  of  corrosion  damage  and  gap  at  the  second  layer  of  an  aircraft  structure.  Detection 
capabilities  for  slots  0.5  inches  long  have  been  demonstrated  by  experiment  and  by  modeling  for 
an  outside  skin  thickness  of  0.160  inch  (2024  T3  aluminum).  Commercial  aircraft  skin  thickness 
is  nominally  0.037  inch  thus  a  favorable  signal  to  noise  ratio  has  been  demonstrated  for  practical 
application.  Work  remains  to  optimize  /balance  the  detection  margin  as  a  function  of  scanning 
speed.  Two  prototype  design  iterations  of  the  system  have  been  built  and  are  under  field 
evaluation  at  Boeing.  These  iterations  will  be  invaluable  in  addressing  application  conditions 
such  as  rivets,  substructure  and  tear  straps.  In  addition,  the  capability  for  measurement  of 
corrosion  and  gap  is  expected  to  be  a  benefit  to  the  development  of  acceptance  criteria  and 
corrosion  analysis  /  management  tools.  Further  engineering  development  support  are  necessary  to 
produce  a  mature  design  that  can  be  considered  for  licensing,  instrument  production  and 
marketing.  It  is  expected  that  the  instrument  design  and  supporting  technology  may  reach 
sufficient  maturity  in  the  next  year  to  warrant  reconsideration  for  a  Phase  II STTR  program  and 
for  licensing,  commercial  instrument  development  and  marketing. 

The  tangible  outputs  of  this  program  are  the  measurement  modeling  software  that  may  be 
marketed  by  the  Center  for  Nondestructive  Evaluation  and  the  Probability  of  Detection  Macro’s 
that  are  marketed  by  D&W  Enterprises,  LTD.  Both  programs  are  Copyright  ©  by  the  developers 
and  marketing  options  will  remain  with  the  developers. 

Pulsed  eddy  current  technology  has  been  demonstrated  to  have  potential  for  detection  and 
measurement  of  corrosion  in  the  second  layer  of  aircraft  structures.  Significant  progress  in 
understanding  the  technology  and  in  addressing  application  conditions  has  been  made  under  the 
current  project  and  in  conjunction  with  related,  concurrent  project  funding  and  support.  Work 
performed  under  the  current  program  is  considered  to  be  a  success  story  but  did  not  progress  to 
licensing  and  commercial  marketing  stage  during  this  Phase.  It  is  anticipated  that  a  Phase  II 
STTR  program  will  be  pursued  pending  further  development  to  bring  it  to  a  level  for 
consideration  of  license,  instrument  development  and  commercial  production.  It  is  clear  that  the 
need  for  acceptance  criteria,  engineering  analysis  tools  and  establishment  of  requirements  for 
measurement  in  maintenance  operations  are  required  before  commercial  marketing  can  be 
considered.  These  developments  are  inter-related,  but  independent  of  consideration  for  further 
development  of  the  pulsed  eddy  current  technology. 

The  authors  appreciate  the  support  of  the  Air  Force  Office  of  Scientific  Research  and  the  data 
shared  by  the  concurrent  industry  and  FAA  programs  that  have  supported  pulsed  eddy  current 
technology  development.  The  combined  efforts  and  resources  were  of  significant  benefit  to  this 
program. 
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Preface 


This  research  reported  here  is  motivated  by  the  need  to  keep  an 
aging  fleet  of  military  aircraft  at  a  safe  state  of  readiness  with  lim¬ 
ited  capital  resources.  A  key  requirement  is  for  the  development 
of  flexible  inspection  technologies  that  can  be  adapted  to  deal  with 
a  wide  range  of  operational  problems  on  a  variety  of  aircraft  each 
of  which  have  their  own  failure  mechanisms.  Typically  such  prob¬ 
lems  are  unforeseen  and  require  new  inspections  to  be  devised  and 
implemented  at  short  notice.  This  program  seeks  a  cost  effect  way 
of  designing  new  inspections  in  response  to  operational  needs  and 
to  validate  the  reliability  of  those  inspections  in  order  that  aircraft 
are  maintained  in  a  safe  condition  for  their  full  or  extended  service 
life. 

The  need  for  cost-effective  assessment  of  aircraft  inspection  reliabil¬ 
ity  is  being  met  in  this  program  by  the  development  of  computer- 
based  tools  for  the  evaluation  of  flaw  detection  probabilities.  In 
addressing  the  need,  we  have  focused  on  the  reliability  of  pulsed 
eddy-current  methods  for  the  detection  of  corrosion  in  aircraft.  Us¬ 
ing  a  new  approach  for  reliability  assessment,  a  measurement  model 
simulates  the  test  conditions  and  sensor  flaw  interactions,  replacing 
a  costly  and  lengthy  process  of  sample  preparation  and  laboratory 
testing.  The  approach  will  allow  the  performance  of  inspections  to 
be  designed,  compared  and  optimized  using  a  technology  that  is 
flexible,  rapid  and  economic. 
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Chapter  1 


Requirements,  Aims  and  Objectives 


1.1  Aims 

The  general  long  term  aim  of  the  present  program,  guiding  the  work  in  both  Phase  I  and 
Phase  II,  is  to  develop  quantitative  inspection  methods  that  are  integrated  with  a  set  of 
tools  and  procedures  to  provide  performance  guarantees  for  first  time  users  and  permit  new 
inspections  to  be  evaluated  in  terms  of  their  performance  reliability.  Inspection  reliability 
is  quantified  using  estimates  of  the  probability  of  detection  (POD)  and  is  dependent  on 
defect  parameters  and  measurement  conditions. 

The  work  program  has  been  subdivided  into  tasks  design  to  achieve  a  number  of 
specific,  well  defined  objectives  in  support  of  the  overall  program  aim.  Tasks  and  Phase 
I  objectives  are  discussed  in  the  next  section. 

The  tools  developed  to  assess  inspection  performance  and  reliability  are  geared  specif¬ 
ically  to  an  emerging  technology;  pulsed  eddy-current  testing.  A  pulsed  eddy-current 
inspection  system  has  been  developed  in  an  earlier  AFOSR  program  at  the  Center  for 
Non-destructive  Evaluation,  Iowa  State  University.  With  additional  support  from  the 
Federal  Aviation  Administration,  the  technology  has  reached  the  point  at  which  beta  test 
equipment  is  being  designed  and  constructed  for  delivery  to  users  in  the  civil  sector.  In 
Phase  II  of  the  current  program,  an  advanced  inspection  system  will  be  designed  and  built 
for  application  to  military  aircraft.  In  Phase  I,  the  focus  has  been  on  the  development 
of  software  to  assess  performance  reliability  and  to  simulate  measurements  of  corrosion 
damage. 

The  central  aim  of  Phase  I  is  to  develop  a  measurement  model,  validated  using  exper¬ 
iment  data  and  coupled  with  a  basic  POD  evaluation  capability  to  provide  a  reliability 
metric  for  corrosion  damage  assessment  on  aircraft.  A  schematic  representation  of  data 
and  model  priorities  is  shown  in  Figure  1.  Inputs  to  the  measurement  model  are  flaw 
and  inspection  parameters  including  physical  parameters  defining  the  sensor  properties, 
the  flaw  size,  shape  and  location  together  with  relevant  material  parameters.  Outputs 
are  the  flaw  dependent  probe  signal  response  as  a  function  of  probe  position,  frequency 
in  the  case  of  a  time-harmonic  excitation  or  time  in  the  case  of  a  transient  response.  The 
calculated  signal  is  used,  together  with  empirical  data  on  the  variability  of  the  inspection 
parameters,  to  estimate  the  probability  of  detection  of  a  given  flaw  and  thereby  assess  the 
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Figure  1.1:  Schematic  diagram  showing  the  measurement  and  POD  models  together  with 
data  input/output  files  and  documents. 

inspection  performance. 

Although  the  value  of  reliability  metric  for  management  of  nondestructive  inspection 
(NDI)  is  widely  recognised,  traditional  methods  of  establishing  detection  probabilities 
through  multiple  sample  preparation  and  extensive  testing  can  be  prohibitively  expensive 
and/or  impractical.  Substitution  of  a  measurement  model  is  a  cost-effective  rapid-response 
alternate  to  the  traditional  approach  [1],  Clearly,  it  is  vital  that  evaluations  of  sensor 
output,  here  performed  through  analytical  and  numerical  treatment  of  the  field  equations 
governing  probe  flaw-interaction,  are  accurate.  Ultimately  it  is  necessary  to  simulate 
realistic  part  geometries  taking  account  of  common  structural  features  of  an  aircraft,  such 
as  installed  fasteners,  lap  joints  and  spars.  However,  initial  progress  has  been  made  by 
analyzing  simply  geometrical  forms  for  the  test  pieces.  In  Phase  II,  the  measurement 
model  will  be  generalize  to  deal  with  complex  structures. 

1.2  Program  Tasks  and  Objectives 

The  proposal  for  the  present  program  identified  three  Phase  I  Tasks  in  line  with  the  need 
to  develop  quantitative  NDE  capabilities  for  life-cycle  management  of  corrosion  in  aircraft 
structures: 

1.  Development  of  quantified  flaw  detection  requirements  and  acceptance  criteria. 
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2.  Development  of  measurement  modeling  tools  to  enable  prediction  of  output  signal 
levels  with  different  probe  designs  and  inspection  conditions. 

3.  Development  and  linking  of  measurement  output  to  a  predictive  POD  model  to 
support  NDE  engineering  especially  in  facilitating  the  introduction  and  validation 
of  new  procedures. 

Here  the  outcomes  of  Task  2,  the  measurement  model  for  performance  and  reliability 
assessment  is  described. 

Measurement  models  can  be  constructed  using  a  number  of  different  approaches  which 
can  be  classified  as  synthetic  or  deterministic.  In  the  synthetic  approach  the  model  is 
a  secondary  data  source  acting  to  integrate,  organize  and  process  data  delivered  to  it 
from  a  primary  source.  Statistical  techniques  or  artificial  neural  networks  fall  into  this 
category  because  they  require  extensive  input  of  empirical  data  in  order  to  be  able  to 
regenerate  probe  signals  in  agreement  with  laboratory  experiments  or  field  inspections.  It 
has  already  been  emphasized  that  reliance  on  sample  preparation  and  laboratory  testing 
for  controlled  sensor  data  is  expensive  and  time  consuming,  severely  limiting  the  ability 
to  quantify  reliability.  Statistical  or  artificial  neural  network  methods  relying  solely  on 
empirical  input  data  suffer  from  the  same  resource  limitations.  Furthermore  they  are 
lacking  a  genuine  predictive  capability  and  therefore  cannot  be  used  beyond  the  limits  of 
the  primary  data. 

Methods  that  ignore  the  physical  relationship  between  probe  and  defect  can  be  pa¬ 
rameterized  or  trained  to  conform  with  the  physical  realities  using  a  deterministic  model 
but  then  the  latter  is  a  prerequisite  and  takes  priority.  However,  the  availability  of  a 
deterministic  model  does  not  render  the  synthetic  methods  redundant  since  these  may 
have  advantages  such  as  speed  of  response  and  ease  of  use. 

Eddy-current  signal  calculations  have  been  carried  out  in  this  program  using  physical 
(deterministic)  models  based  on  electromagnetic  field  theory.  Limited  testing  of  the  output 
via  comparisons  with  experiment  are  necessary  but  only  a  few  experiments  are  needed 
for  validation  purposes  and  this  process  is  therefore  not  expensive.  This  approach  is  the 
only  effective  solution  that  avoids  the  cost  of  multiple  sample  preparation  and  extensive 
empirical  testing.  A  new  measurement  model  based  on  the  physics  of  the  probe-defect 
relationship  has  been  developed  for  Task  2.  The  theoretical  study  of  the  interaction  of 
eddy  current  probe  and  flaw  yields  a  set  of  equations  that  are  solved  numerically  to  give  the 
electromagnetic  field  in  the  presence  of  corrosion  damage,  the  corrosion  being  represented 
as  a  non-conducting  region  embedded  in  a  conducting  region.  The  field  equations  have 
been  translated  into  an  efficient  computer  code  using  numerical  methods  and  solutions 
generated  from  input  data  which  describe  the  test  conditions. 

The  theory  of  transient  electromagnetic  fields  has  been  used  to  evaluate  pulsed  eddy- 
currents  (PEC)  in  conductors.  By  computing  numerical  solutions  of  the  field  equations, 
probe  signals  due  to  localized  regions  of  corrosion  have  be  determined  for  the  first  time  in 
this  study  as  a  function  of  probe  position  and  time.  Previous  theoretical  studies  of  pulsed 
eddy-current  nondestructive  evaluation  have  represented  corrosion  as  an  infinite  layer  [2, 
3,  5].  In  another  earlier  study,  pulsed  eddy-current  interactions  with  cracks  have  been 
predicted  and  compared  with  experimental  measurements  [4],  The  new  results  presented 
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here,  however,  are  derived  from  a  more  flexible  model  than  has  been  used  previously.  The 
model  makes  use  of  a  regular  three-dimensional  array  of  volume  elements  which  can  be 
assembled  in  a  form  that  represents  either  a  crack  or  a  finite  region  of  corrosion.  Simply  by 
assigning  conductivity  values  to  the  volume  cells  in  a  three  dimensional  grid,  a  corroded 
region  arbitrary  shape  can  be  represented,  including  for  example  an  approximation  of  a 
conical  corrosion  pit. 

Task  2  entails  the  following  objectives  which  have  been  accomplished  with  the  aid  of 
transient  eddy-current  measurements  models: 

1.  An  existing  transient  eddy  current  code  has  been  modified  for  predicting  signals 
due  to  corrosion.  Parametric  studies  of  probes  have  been  carried  out  find  probes 
with  good  signal  sensitivity.  The  corrosion  for  this  study  was  modeled  as  a  uniform 
infinite  non-conducting  layer. 

2.  Noise  and  signal  variability  parameters  have  been  measured  for  a  pulsed  eddy- 
current  test  instrument  and  scanning  system. 

3.  A  response  model  has  been  developed  to  determine  signals  due  to  finite  regions  of 
corrosion. 

4.  Interface  the  measurement  model  to  the  POD  code. 

5.  This  report  has  been  produce. 

The  results  from  Task  2  are  reported  in  next  chapter.  To  conclude  this  chapter,  the 
prospects  for  finding  quantitative  measures  of  corrosion  damage  are  reviewed  below. 


1.3  Quantitative  Corrosion  Assessment 

The  rough  location  of  corrosion  damage  beneath  an  aircraft  skin  can  be  identified  directly 
from  raw  eddy-current  data,  using  for  example,  false  color  c-scan  images.  The  data 
needed  to  generate  such  images  is  acquired  by  scanning  the  surface  of  the  aircraft  using 
an  eddy-current  probe  while  recording  transient  or  time-harmonic  field  measurements.  It 
is  relatively  easy  to  get  an  approximate  idea  of  the  lateral  position  of  the  corrosion  from 
the  images  but  there  is  a  need  to  estimate  material  loss  and  to  determine  the  depth  of 
corrosion  damage  in  a  layered  structure.  This  information  is  not  readily  available  from 
the  c-scan  images.  The  task  of  locating  the  in-depth  limits  of  localized  corrosion  is  an 
important  element  in  the  planned  Phase  II  program  and  can  be  accomplished  with  the 
aid  of  the  measurement  model  developed  in  Phase  I. 

The  simplest  model  of  corrosion  treats  the  test  piece  as  a  series  of  uniform  layers 
stacked  on  top  of  one  another.  The  parameters  of  the  structure,  the  conductivity  and 
thickness  of  the  layers,  has  been  found  from  pulsed  eddy-current  data  by  parameter  fit¬ 
ting  [5].  However,  the  infinite  layer  model  does  not  give  realistic  representation  of  a 
localize  and  irregular  region  of  corrosion  damage.  For  the  case  of  crack  inversion,  it 
has  been  shown  that  irregular  shapes  can  be  reconstructed  using  eddy  current  data  by 
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Figure  1.2:  Schematic  diagram  showing  the  measurement  model  in  an  iterative  inversion 
process. 

applying  a  least  square  optimization  scheme  that  relies  on  a  fast  and  accurate  forward 
measurement  model  to  predict  probe  signals  [6,  7].  By  adjusting  flaw  parameters  itera¬ 
tively  until  calculated  signals  agree  with  measurements,  Figure  1.2,  the  boundaries  of  the 
crack  was  found.  Therefore  we  can  be  reasonably  confident  that  a  similar  approach  would 
determine  the  limits  of  the  corroded  region  from  eddy-current  probe  data.  The  general 
principles  for  accomplishing  this  task  have  been  established  by  Norton  and  Bowler  for 
time-harmonic  eddy-current  data  [8]  but  it  is  not  difficult  to  extend  these  principles  to 
deal  with  transient  fields. 

Implementation  of  the  inversion  procedure  for  the  quantification  of  material  loss  due 
to  corrosion  using  pulsed  eddy-current  data  will  be  set  as  a  Phase  II  objective.  This 
objective  will  be  accomplished  with  the  aid  of  a  measurement  model  that  can  calculate 
pulsed  eddy-current  signals  due  to  local  regions  of  material  loss  and  has  been  developed 
in  Phase  I. 
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Chapter  2 


Simulation  of  Pulsed  Eddy-Current 
Probe  Response  due  to  Corrosion 

2.1  Introduction 

In  order  to  quantify  material  loss  due  to  corrosion  using  eddy-current  techniques,  advances 
are  needed  in  inspection  performance  and  data  analysis.  Fortunately  these  advances  axe 
emerging  with  the  use  of  new  field  sensors  and  measurement  methods  such  as  pulsed 
excitation  or  remote  field  methods.  Laboratory  tests  conducted  at  the  Center  for  Non¬ 
destructive  Evaluation  at  Iowa  State  University  show  that  a  2%  material  loss  could  be 
reliably  detected  under  laboratory  conditions  using  pulsed  eddy-currents.  For  conventional 
eddy-current  techniques,  by  comparison,  the  detection  limit  is  more  than  10%  material 
loss  due  to  simulated  corrosion  in  lap  joints  at  a  probability  of  detection  of  level  90%  [10]. 

In  transient  eddy-current  inspection,  a  current  pulse  in  a  coil  excites  a  time- varying 
electromagnetic  field  which  induces  current  in  a  conducting  test  piece.  As  this  pulse 
migrates  into  the  material  it  is  broadened  by  dispersion  and  perturbed  by  discontinuities 
in  the  conductivity  and  permeability  of  the  conductor.  Sub-surface  defects  disrupt  the 
field  causing  a  perturbation  to  migration  back  to  the  surface.  Above  the  conductor,  it 
can  be  observed  as  a  transient  signal  in  the  time  domain,  either  as  an  electromotive  force 
(emf)  in  a  pick-up  coil,  or  as  a  signal  produce  by  a  solid  state  field  sensor  such  as  a 
Hall  probe  or  magnetoresistive  sensor.  A  clear  transient  signal,  measured  with  the  probe 
remote  from  flaws  or  structural  features  such  as  fasteners,  is  typically  used  as  a  reference 
which  is  subtracted  from  all  other  signals.  The  difference  signal  contains  flaw  information, 
the  fastener  signal  or  indeed  any  indication  which  differs  from  the  reference  signal.  The 
transient  is  broad  band,  carrying  information  from  a  wide  frequency  spectrum  in  contrast 
to  a  time-harmonic  excitation  which  yields  information  only  at  a  single  frequency. 

Eddy-current  inspection  is  effective  for  detecting  surface  or  near  surface  cracks  and 
corrosion  in  conductors.  Because  induced  currents  in  a  plane  structure  flow  roughly 
parallel  to  the  interfaces,  cracks  that  are  predominantly  perpendicular  to  the  surface  are 
readily  detectable,  the  response  being  related  to  the  area  of  the  crack  face.  Corrosion 
in  aircraft  structures,  such  as  exfoliation  corrosion  and  corrosion  in  lap  joints,  tends  to 
emerge  at  the  interfaces  between  the  plates  that  make  up  the  structure.  Eddy  current 
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Figure  2.1:  Hall  sensor  mounted  on  the  axis  of  an  excitation  coil. 


measurements  are  affected  by  the  change  in  conductivity  caused  by  decomposition  of  the 
metal  but  the  presence  of  corrosion  is  not  easy  to  detect  at  the  early  stages  due  to  the 
small  perturbations  of  the  induced  current. 


2.2  Measurement  Model  of  Pulsed  Eddy- Currents 

In  order  to  assess  the  capabilities  of  eddy  currents  for  detecting  and  characterizing  cor¬ 
rosion,  the  uniform  layer  model  has  been  used  as  an  idealized  representation  of  corrosion 
attack  between  plates.  The  probe  signal  is  initially  calculated  for  the  frequency  domain 
using  analytical  formulae  and  then  transformed  to  the  time  domain  using  a  fast  Fourier 
transform. 

In  experiments  on  a  simple  layer  structure,  a  series  of  electromagnetic  field  pulses 
were  induced  in  the  specimen  by  passing  a  square  wave  current  through  an  excitation 
coil.  Transient  fields  interact  with  the  material  and  with  the  non-conducting  regions 
between  plates  giving  rise  to  a  field  perturbation  which  migrates  to  the  surface.  The  near 
surface  field  is  measured  as  a  function  of  time  using  a  Hall  sensor,  Fig.  2.1.  The  pulsed 
eddy-current  layer  model  has  been  evaluated  by  comparing  predictions  with  experiment, 
Figure  2.2.  The  results  show  that  three  distinct  transient  signals  axe  found  depending  on 
whether  the  air  gap  between  two  plates  is  due  to  loss  of  material  from  the  top  plate,  loss 
of  material  from  the  bottom  plate  or  simply  an  increase  in  plate  separation [2]. 

The  layer  model  is  adequate  for  basic  probe  optimization  studies  and  for  finding  an 
upper  bound  to  the  probability  of  detection  of  a  region  of  corrosion.  However  it  is  limited 
because  it  cannot  be  used  to  compute  the  effects  of  localized  corrosion  damage.  In  order 
to  perform  such  calculations  for  finite  pockets  of  material  loss,  a  new  volume  element  code 
has  been  developed.  This  calculates  probe  signals  directly  in  the  time  domain  making  use 
of  an  integral  equation  with  a  time-domain  kernel  [9].  The  volume  element  approach  uses 
the  same  kernel  as  one  previously  used  for  boundary  element  calculations  and  produced 
the  same  results  when  used  to  simulate  the  effects  of  a  subsurface  crack,  Fig.  2.3. 


Plate  thickness  2.00mm,  gap  0.45mm. 


Figure  2.2:  Comparison  of  prediction  of  a  layer  model  of  corrosion  with  experiment. 

2.3  Phase  I  Accomplishments 

The  accomplishments  below  correspond  to  the  objectives  given  in  Section  1.2. 

2.3.1  Investigation  of  Probe  Sensitivity 

The  sensitivity  of  probes  have  been  compared  using  on  the  basis  of  signal  to  noise  ratio.  A 
probe  containing  a  Hall  sensor  gives  a  magnetic  field  measurement,  H ,  (in  A/m)  normal¬ 
ized  to  the  excitation  coil  current.  In  general,  the  axial  length  and  the  inner  radius  of  the 
coil  is  small  for  optimum  sensitivity.  The  coil  lift-off  (height  of  the  base  above  the  work- 
piece)  should  be  as  small  as  possible  and  the  Hall  sensor  placed  as  close  to  the  surface 
of  the  work-piece  as  possible  on  the  coil  axis.  The  sensor  is  supplied  in  an  encapsulated 
dual-in-line  package  whose  size  constrains  the  probe  and  coil  lift-off.  It  was  concluded 
that  the  probe  dimensions  for  best  performance  are  those  given  in  Table  2.1. 

2.3.2  Noise  and  signal  variability 

Quasi-random  background  noise  levels  for  the  above  probe  was  approximately  25  A/m  per 
Ampere  of  coil  current  whereas  signal  variability  (standard  deviation)  was  approximately 
100  A/m,  mainly  due  to  variations  of  lift-off  (i.e.  inconsistency  in  the  probe-workpiece 
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Table  2.1:  Coil  Parameters. 


Outer  radius 
Inner  radius 
Axial  length 
Nominal  lift-off 
Number  of  turns 


10.0  ±  0.05  mm 
2.50  ±  0.01  mm 
2.80  ±  0.01  mm 
0.80  ±  0.01  mm 
1000  ±  1  mm 


contact).  The  signal  variation  compares  with  the  typical  flaw  signal  levels  shown  in  figure 
2.2  for  a  material  loss  of  22.5%.  Note  that  a  1%  material  loss  corresponds  to  a  signal 
level  of  80  A/m  which  is  about  three  times  the  noise  level.  This  is  roughly  the  limit  of 
resolution  of  material  loss. 

Signal  variability  or  lack  of  repeatability  is  mainly  dependent  on  scanning  conditions 
especially  the  positioning  of  the  probe  on  the  surface  of  the  work-piece.  The  main  cause 
of  the  lack  of  repeatability  in  the  signal  is  probe  lift-off  or  probe  tilt.  Assuming  the  signal 
variability  is  due  to  a  lift-off  only,  it  is  possible  to  calculate  the  lift-off  variation  which 
give  rise  to  the  observed  variation  of  the  signal.  In  this  way  the  lift-off  variation  of  0.15 
mm  was  inferred. 


2.3.3  Response  Model  for  Localized  Regions  of  Corrosion 

The  measurement  model  for  calculating  the  probe  response  due  to  local  pockets  of  corro¬ 
sion  of  arbitrary  shape  is  described  in  Chapter  3. 

2.3.4  Calculation  of  Probability  of  detection 

A  flaw  is  said  to  be  detected  or  detectable  when  the  signal  it  produces  v,  exceeds  a  given 
threshold,  vt.  The  threshold  level  is  chosen  as  a  compromise  between  maximizing  the 
probability  of  detection,  which  means  vt  should  be  low,  and  minimizing  the  number  of 
false  calls  which  requires  a  high  value.  The  threshold  can  be  calculated  from  a  statistical 
analysis  of  the  observation  data  based  on,  for  example,  a  tolerable  level  of  false  calls,  but 
for  present  purposes  we  shall  simple  set  the  detection  threshold  at  three  times  the  noise 
level.  This  seems  to  be  a  satisfactory  criterion  for  most  purposes. 

Repeated  measurement  produces  a  distribution  of  values  of  the  signal  v  about  a  mean 
value  say  v(x),  with  a  standard  deviation  of  av.  Assuming  this  distribution  is  normal 
(Gaussian),  the  probability  that  a  given  signal  for  a  flaw  of  size  x  exceeds  the  threshold 
and  therefore  satisfies  the  criterion  for  detection,  is  given  by 

1  roc 

x)  =  nr  I  exp[-(t>-e(aO)  /2av(x)]dv 
y J2irav  Jvt 

=  ^erfc(V)  (2.1) 

where  V  =  (vt  —  v{x))/av  and  erfc  is  the  complementary  error  function. 

It  is  reasonable  to  suppose  that  the  mean  value  of  the  observation  is  approximated  by 
the  model  prediction  for  a  given  flaw  of  characteristic  dimension  x.  Applying  the  layer 
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model  to  the  calculation  of  POD,  with  x  representing  the  reduction  in  plate  thickness  due 
to  loss  of  material,  the  corresponding  detection  probabilities  are  shown  in  Fig.  2.4  for 
cases  where  the  first  layer  plate  is  2  mm  and  4  mm  respectively.  . 


2.4  Volume  Element  Measurement  Model  Descrip¬ 
tion 

The  model  input  is  a  set  of  inspection  parameters  used  to  characterize  the  physical  condi¬ 
tions  of  an  idealize  inspection.  These  parameters,  which  include  the  probe  dimensions,  the 
material  properties  of  the  part  and  the  flaw  dimensions,  are  used  to  calculate  the  sensor 
probe  response  to  a  flaw,  such  as  a  region  of  corrosion.  The  flaw  signals,  the  variability 
of  the  input  parameters,  noise  sources  and  the  detection  threshold  are  to  be  combined  in 
a  probability  of  detection  model  to  determine  inspection  performance. 

An  eddy-current  probe  responds  to  the  loss  of  conducting  material  due  to  corrosion 
but  it  also  detects  probe  lift-off  changes,  fasteners  and  other  structural  features  of  the 
aircraft.  In  order  to  take  account  of  these  effects,  the  measurement  model  has  the  capa¬ 
bility  of  representing  general  three-dimensional  objects.  However,  the  time  taken  for  the 
calculations  is  kept  within  reasonable  bounds.  In  developing  a  modern  three-dimensional 
electromagnetic  field  code,  much  of  the  software  engineering  effort  is  devoted  to  reducing 
the  computation  cost  of  the  calculation. 

Initially,  the  model  has  been  used  to  represent  a  finite  pocket  of  corrosion  in  an  other¬ 
wise  uniform  plate.  However  the  approach  taken  is  quite  general  and  can  be  adapted  to 
accommodate  fasteners  and  lap  joints  etc.  without  excessive  computational  cost.  These 
features  will  be  included  as  part  of  the  Phase  II  extensions.  The  probe  signal  is  evaluated 
from  the  electromagnetic  field  in  the  part,  the  field  being  calculated  using  a  new  volume 
element  method  based  on  recent  advances  in  time-domain  electromagnetic  theory  [9].  A 
region  of  corrosion  is  represented  as  a  non-conducting  or  partially  conducting  zone  of 
arbitrary  shape  constructed  using  volume  elements. 

The  measurement  model  code  determines  the  eddy-current  field  in  the  part  using  a 
matrix  solver.  This  is  a  three-stage  calculation  in  which  matrix  elements  are  calculated  in 
the  first  stage  and  in  the  second,  the  field  of  the  probe  is  determined.  In  the  third  stage, 
the  field  in  the  flaw  is  found  by  solving  a  matrix  equation  and  the  probe  signal  calculated 
from  the  solution.  Details  of  the  field  analysis  and  numerical  method  are  given  in  Chapter 
3. 
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Figure  2.3:  Calculations  of  the  response  of  a  field  sensor  on  the  axis  of  a  coil  in  the 
presence  of  a  sub-surface  crack  or  notch.  The  vertical  scale  gives  the  magnetic  field  at  the 
site  of  the  Hall  sensor  (A/m).  The  double  peak  corresponds  to  the  region  where  the  notch 
is  approximately  1  mm  from  the  surface.  To  the  left  of  the  double  peak  is  the  indication 
from  the  deeper  portion  of  the  notch.  This  part  of  the  indication  is  delayed  in  time. 


Chapter  3 

Numerical  Evaluation  of  Pulsed 
Eddy-Current  Interaction  with 
Volumetric  Flaws 


In  the  problem  considered,  induce  current  pulses  excited  by  an  eddy-current  probe  in¬ 
teract  with  a  three  dimensional  flaw  in  an  otherwise  uniform  half-space  conductor.  The 
perturbation  of  the  electromagnetic  field  due  to  the  flaw  is  detected  by  measuring  the  time 
variation  of  the  external  magnetic  field  or  the  induced  emf  in  a  sensor  coil.  The  aim  is  to 
compute  the  field  as  a  function  of  time  and  hence  predict  the  sensor  signals.  A  volume 
integral  formulation  of  the  problem  has  been  developed  in  order  to  calculate  the  perturbed 
field  due  to  a  region  of  corrosion  or  a  crack  in  an  otherwise  homogeneous  conductor.  An 
integral  equation  for  the  field  in  the  flaw  has  been  derived.  It  is  then  approximated  by  a 
matrix  equation  using  the  moment  method.  Closed  form  expressions  have  been  obtained 
for  calculating  the  matrix  elements  including  the  singular  element.  From  the  solution  of 
the  matrix  equation,  the  external  magnetic  field  can  be  evaluated  using  a  formula  derived 
using  reciprocity  principles. 

3.1  Introduction 

The  evolution  of  the  electromagnetic  field  in  a  conductor  containing  an  inhomogeneity 
has  been  calculated  using  volume  elements  by  San  Filipo  and  Hohmann[12]  and  by  Wang 
and  Hohmann[13]  using  the  finite  difference  technique.  These  calculations  are  aimed  at 
applications  in  geophysics  but  similar  problems  exist  in  pulsed  eddy-current  nondestruc¬ 
tive  evaluation (NDE).  In  NDE,  the  inhomogeneity  is  usually  considered  as  localised  which 
favours  the  volume  element  technique  because  its  primary  goal  is  calculate  the  field  in  the 
flaw  region  only.  If  this  region  is  small  and  well  defined  then  a  good  approximation  may 
be  obtained  using  a  relatively  small  number  of  volume  elements. 

This  article  describes  a  time-domain  volume  element  formulation  for  obtaining  solu¬ 
tions  that  are  approximated  as  piecewise  linear  in  time  and  space.  Results  of  transient 
field  calculations  are  not  included  at  this  stage.  However,  the  key  features  of  the  approach 
are  set  out,  including  the  calculation  of  the  self-element  taking  into  account  the  singularity 
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of  the  Green’s  function.  In  addition,  the  remaining  matrix  elements  corresponding,  to  a 
free-space  dyadic  Green  kernel  have  been  expressed  in  closed  form.  From  the  information 
given,  it  should  be  a  relatively  straightforward  matter  to  adapted  the  matrix  calculation 
to  other  geometries  such  as  a  conducting  half-space. 


3.2  Formulation 


Assuming  that  the  material  properties  are  linear,  that  the  permeability  is  that  of  free 
space  and  that  displacement  current  is  negligible,  the  electromagnetic  field  in  a  conductor 


containing  a  local  conductivity  variation  satisfies 

v  V  EM  -  ,/ 

V  x  E(r)  =  Ho  dt 

(3.1) 

and 

V  x  H(r)  =  (T0E(r)  +  P(r) 

(3.2) 

where  a0  is  the  host  conductivity  and  P(r)  =  [a(r)-a0]E(r)  is  a  current  dipole  distribution 
representing  the  effect  of  the  flaw  as  an  equivalent  source. 

A  solution  of  Maxwell’s  equations  will  be  sought  using  the  magnetic  vector  potential 
by  following  standard  electromagnetic  theory.  The  vector  potential  is  defined  such  that 
B  =  V  x  A,  hence,  for  a  medium  with  the  permeability  of  free  space 


H  =  —V  x  A.  (3.3) 

Ho 

By  substituting  for  /ioH(r)  in  the  induction  law,  one  concludes  that  the  electric  field  has 
the  form 

E=-^-v^  (3-4) 

where  Vip  is  a  scalar  potential.  Substituting  (3.3)  and  (3.4)  into  (3.2)  and  choosing 
the  gauge  where  ip  =  — V  •  A//i0  it  is  found  that  the  vector  potential  satisfies  a  three 
dimensional  vector  diffusion  equation 


dA 

V2A  +  Wo-Qj:  =  -to  P(r)- 


(3.5) 


E 


=  vv)  -  a, 

V  at  lio<j0  ) 


(3.6) 


Combining  the  gauge  condition  with  (3.4),  shows  that  the  electric  field  is  given  by 

d_ 

dt  nQcr0 

where  1  =  xx  +  yy  +  zz  is  the  unit  dyad.  The  solution  of  (3.5)  will  be  written  in  the  form 

A(r)  =  A(0)(r)  ~  Ho  Jq  ^^(r,r',A  -  t1)  ■  P(r \£)  dr1  dt1.  (3.7) 

where  A^°l(r)  is  the  unperturbed  field  and  the  dyadic  Green’s  function  Qa{t,  r',  t  —  t1) 
satisfies  ; 

V2^(r,r',t  -  t’)  +  l*oa0^^1^- — —  =  -15 (r  -  r')6(t  -  t').  (3.8) 

A  solution  of  (3.8)  for  a  half  space  conductor  are  given  in  reference  [9]. 
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3.3  Electric  field  integral  equation 


From  (3.6)  and  (3.7),  the  electric  field  in  the  presence  of  a  volumetric  flaw  is  written  as 

E(r,  t)  =  E<°> (r,  t)  -  _  i-vv)  •  jf*  Jqa{ r,  r',  f,  t')  •  P(r',  t')  dr'  dt',  (3.9) 

where  P(r',  t')  is  the  induced  dipole  density  at  the  defect  and  is  the  flaw  region.  Mul¬ 
tiplying  by  a(r)  -a0  =  a0v(r)  gives 

P(r,  t)  =  P<°)  (r,  t)  -  a0v(r)  w)  •  j*  JQa(t,  r',  t,  t)  •  P(r',  H)  dr'  dt',  (3.10) 

which  is  a  Fredholm  integral  equation  of  the  second  kind  for  the  flaw  polarisation.  In 
Equation  (3.10),  P(r ,t)  =  aoi>(r)E(r,£)  and  P(0)(r,t)  =  <70^(r)E^(r,  f). 


3.4  Integral  kernel 

In  the  case  of  an  unbounded  domain,  the  integral  kernel  is  given  by [9] 

Qa{  r,  =l(j)(\r-r'\,t-t'), 


(3.11) 


where  <f>(\r  —  r'\,t  —  f)  satisfies  the  scalar  diffusion  equation  with  a  singular  source  term, 
the  solution  being 

m,t  - 1)  =  -  t).  (3.12) 

where  H(t  —  t')  is  the  Heavyside  step  function,  unity  for  positive  argument  and  zero 
otherwise.  This  scalar  Green’s  function  can  also  be  written  as 


a(t  -  t ')  =  2 \fiT-  t')/fioaQ 


where 

is  a  diffusion  length  and 


U(s)  =  ~e~‘\ 


Note  that  the  function  U  ( s )  is  normalised  in  the  sense  that 

/OO 

U(s)  ds  —  1. 

-OO 

An  alternative  to  equation  (3.12)  is  to  write 

/W(|r  -  r'|,i  -  ()  =  Ajtl(£)ff(t  -  f). 


(3.13) 

(3.14) 

(3.15) 

(3.16) 

(3.17) 
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The  integral  of  the  scalar  Green’s  kernel  is  of  interest,  hence  we  define 

=  (3-18) 

Note  that  as  t  —  t'  -*  oo,  the  complimentary  error  function  approach  unity  and  the  right 
hand  side  of  (3.18)  tends  to  the  familiar  static  Green’s  function. 

In  order  to  make  prediction  that  can  be  compared  with  experiment,  it  would  be  nec¬ 
essary  to  used  a  half-  space [9]  or  plate  kernel  in  (3.10).  The  addition  complications 
associated  with  this  will  not  be  considered  here.  Instead,  the  procedure  of  discretisation, 
as  prescribed  by  the  moment  method,  is  defined  for  a  flaw  in  an  unbounded  domain. 


3.5  Regularisation 

The  differential  operators  might  be  dealt  with  in  a  numerical  scheme  using  difference  for¬ 
mulae  but  a  better  approach  is  to  bring  the  derivatives  inside  the  integrals  with  respect  to 
time  and  space  allowing  them  to  act  directly  on  the  kernel.  This  leads  to  a  modified  kernel 
that  can  be  approximated  by  adopting  a  numerical  integration  scheme.  Note  however, 
that  the  action  of  the  time  derivative  on  the  time  integral  leads  to  a  form 


(319) 

evaluation  of  which  requires  the  evaluation  of  F(t,  t')  at  t'  =  t.  If  it  happens  that  F(t,  if) 
cannot  be  evaluated  at  t'  =  t,  them  some  other  way  must  be  found  to  bring  the  time 
derivative  inside  the  integral. 

In  the  case  of  the  integral  equation  (3.10),  for  the  dipole  distribution  P(r,  t)  one 
identifies  F(t,  t')  with  the  components  of 


F(r ,t,t')  =  J^gA(r,r',t,t')  •P(r',t')  dr' 


(3.20) 


Assuming  the  dipole  density  is  smooth,  the  integral  can  be  approximated  for  small  t  —  t' 
by 

F ■  *'P(r.(')  (3.21) 

JO-o 

where  is  a  sphere  of  radius  Rq  about  r.  Integration  gives 


F(r,(.f)  =  —  P(r,f)  (3.22) 

Wo 

Provided  P(r,t)  is  smooth,  the  function  F(r,  t,  t')  is  non-singular,  hence  we  can  bring  the 
time  derivative  inside  the  time  integral.  Bringing  it  inside  the  spatial  integral  is  another 
matter  and  does  in  fact  result  in  a  kernel  that  has  a  singularity  in  time  at  t  —  t!  =  0.  In 
finding  a  discrete  form  of  the  equation  that  determines  the  approximate  dipole  density  of 
the  source,  this  singularity  is  dealt  with  by  subtraction.  Thus  we  write 
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where 


P(r,  t)  -  G0v{r)jjZ(r, r',  t)  ■  P(r',  t)  dr'  = 

P(0)M  +  a,v{v)fjg{r,r\t,l!)-[V(T\i!)-V(v\t))dv'dt'  (3.23) 

K(Ty,t)  =  -»0gA(r,r',t)  +  VV-  J gA(\r-r'\,t-t')d(t-t'),  (3.24) 

g(r,r',t-t')  =  -^I-^w)^(r,r't-f')  (3.25) 

For  an  unbounded  domain,  Equation  (3.11)  holds  which  means 

IC( r,r',t)  =  Ko(r,r i)  =  r'|,()I+  -VVSflr  -  r'|,<).  (3.26) 

*0 

and 

g(r,  r',  t-t')  =  g0(  r,  r',  =  -  -  ^w)x0(r,  r',  t  -  t')  (3.27) 

These  functions  are  now  expanded  and  the  integral  operators  of  equation  (3.23)  approxi¬ 
mated  by  a  discrete  form  using  the  moment  method. 


3.6  Expansion  of  the  unbounded  domain  kernel 


In  expanding  equation  (3.27)  it  is  necessary  to  consider  the  time  and  spatial  derivative. 
Firstly  we  have 

dHR,t)- 2 -Kt)’ -*]**■*) 


Mo  O’O' 


dt  a 2 


(3.28) 


Note  that  this  involves  a  t  1  singularity.  Secondly 

d2(j)(R,t)  2  r  /,r,N2 


dp* 


-IKS)  -']«*•* 


Hence 


a0x  •  £/(r,  r',  t-t')-x  =  a0Ggx(r,  r',  t,  t')  =  — 


or 


(v  -  y'Y  +  (z-  z'Y 


a * 


<t>(\ * 


(3.29) 

(3.30) 

V| 

(3.31) 


and 


<r„i  •  e(r.r',(  -  0  ■  y  =  <r0G7(r,r',M')  =  4(x  ^  g'V(|  r  -  r'  |,t  - *')  (3.32) 


a1 


where  x  and  y  are  unit  vectors.  Similar  terms  are  obtain  for  yy  zz  yx  yz  and  zy  components 
by  a  cyclic  rotation  of  the  variables.  The  functions  defined  by  Equation  (3.31)  and  (3.32) 
lead  to  a  t~l  singularity  in  the  time  integral.  Hence  the  need  to  use  some  method  such  as 
the  subtraction  scheme,  to  regularise  the  time  integral. 
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3.7  Primary  Matrix  Evaluation 


The  discrete  form  of  the  equation  for  the  dipole  density  is  found  by  assuming  that  the 
solution  may  be  approximated  as  piecewise  linear  in  space  and  time.  Matrix  elements 
are  then  defined  by  point  matching.  This  process  is  carried  out  in  two  stages.  Firstly 
the  solution  is  approximated  as  piecewise  linear  in  space  by  expanding  in  terms  of  linear 
expansion  functions.  Point  matching  then  produces  a  set  of  time  dependent  primary 
matrix  elements.  Secondly  the  time  variation  is  approximated  as  piecewise  linear. 

Linear  expansion  in  space  and  point  matching  gives  a  set  a  matrix  elements  defined 

by  , 

G<$>(()=/  (3.33) 

where  x,y  =p,  q.  The  function  ^(r ')  which  is  a  convolution  of  the  expansion  and  testing 
functions  is  given  by 

=  'l’j(j-)M-f)M-r)  (3-34) 

0X  Oy  0Z 

where 

Ms)  =  1  ~  \s  ~  j\  (3-35) 

for  j  —  1  <  s  <  j  +  1  and  zero  otherwise. 

For  an  unbounded  domain, 

Qap(t)  =  =  [  Go(ra,r' dr'  =  [  £o(0,r',  t)^-a(r')  dr',  (3.36) 

where  £70(r,  r',  t)  is  given  by  (3.27).  The  components  of  the  tensor,  written  as 

G%(t)=P-e%  («)•«.  (3.37) 

For  the  self  element,  equations  (3.31)  and  (3.36)  combine  to  give 


8  L  Jo*  Jo  a2*  ^)^00(>(r)  dx  dy  dz 

8  fa  rsv  rs*  1  /  y2  +  z2\ 

(t  —  t')  Jo  Jo  Jo  ex3  \  a2  ) 

U(x/a)U(y/a)U(z/a)tpooo{r)  dx  dy  dz 

8  fSx/a  fSy/a  r&z/a 

(TT)L  L  /»  ^-y 

U(x)U(y)U(z)ipooo(ox)  dx  dy  dz 
2  r‘-%_o*\u{x)dxx 

_ — _ p{ 0)  ( — )  f F(2)  ( — )  F(0)  ( — )  + 

SxSySz(t  —  If)  Vcx/L  \a)  \a) 

i?(o)  ^^(2)  _  Fw  ^ijF(o)  ^.y 


(3.38) 


18 


where 


and 


F^n\s)  =  -7=  [  xn(x  —  s)e  x2dx 
y/lX  Jo 

For  the  cases  where  n  =  0, 1, 2,  it  is  found  that 

F^(s)  —  -]=  -  s  erf(s)  —  U(s), 

V71" 

p{l){s )  =  --^=  +  ^erf(s) 

FW(s)  =  4=  ~  Is  erf  (s)  -  U(s). 

V7T  2 

For  a  general  matrix  element,  equations  (3.31)  and  (3.36)  give 

rxx  (i  _  fn  =  [{j+1)Sx  [{k+1)Sv  [{l+1)Sz  _£/\  _  y2  +  z2\ 

J (j - 1)5 x  J(k-l)5y  J(l-l)5z  a2  \  a2  ) 

(j>{ 0,  r,  t  -  t')ipjki(T )  dx  dy  dz 

1  y(j+l)^/Q  i‘(k+i)6v/a  i-(l+l)5z/a ^  2  ^ 

(t  -  t')  J(j-l)5x/a  J(k-l)5v/a  J(l-l)5z/a  ^ 

U(x)U(y)U (z)xpjki(ar)  dx  dy  dz. 


(3.39) 

(3.40) 

(3.41) 

(3.42) 


(3.43) 


This  can  be  expressed  in  terms  of  standard  functions  as  defined  above  by  separating  the 
integrals  as  follows. 

r(k+\)6v/a  r{l+l)5z/a 

J(k- !)*./«  7(J-l)J./a  ^  ^  3 

(1  - 1^  -  *1)  dy dz 


or 


7^(0)  /  °X  N  771(2)  /  °y  \  771(0)  /  °Z  X 

8 SJMt-t')  ‘  'a1!  *  V  ‘  V  + 

t  \  -w-if  O')  /  ^7/  \  -rif  O')  / 


Fk°\-)Fl2]{-)  -  Fk0)(-)Fl{0\-) 

*  v  cr  vcx/  *  K  a J  1  KaJ  J 


a 


(3.44) 


where 


Fjn\s )  =  -7=  f  xn(\x  —  js\  —  s)e  X*dx 
=  (sj+l)  -  2 F^ (Sj)  +  F^ ( Sj _0 

with  Sj  =  js.  Note  that  FQn\s )  =  2F^(s). 


(3.45) 
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Prom  (3.32)  it  can  be  seen  that  the  terms  off  the  tensor  diagonal  are  of  the  type 

fix  rSy  rSz  4 xy 


Goooit-t')  =  8 JQX  jj  JQZ  -^-(p(0,T,t-t')^ooo(T)dxdydz 

ft  rdx/at  rSy/a  rSz/ot 

=  t. — ttt  /  /  /  xyU(x)U(y)U(z)^mo{ar)  dx  dy  dz 

(t  —  v)  Jo  Jo  Jo 


a 


8x5ySz(t-t')  V  V  V 


(3.46) 


For  a  general  off-diagonal  element 

r{j-\-l)Sx  f(k+l)6y  r{l+l)6z  4:XV 

G%(t  - 1')  =  [  f  [  —^<f)(0,T,t-t')ipjki(T)dxdydz 

3  J{j-l)Sx  J(k-l)5y  J{l-l)6z  a* 

1  r(i+ l)Sx/a  Ak+ygy/a  Al+l)6z/a 

=  7 - — ttt  /  /  /  xyU(x)U{y)U(z)xJjjki(r)  dx  dy  dz 

( t  -  t')  J(j-l)Sx/a  J(k-l)6,,/a  J(l-l)Sz/a 


(t  -  t ')  J(j-l)Sx/a  J{k-l)Sv/a  J(l- l)Sz/a 
=  "8 WM)'1  (a)F‘  (a)F'  {a> 

3.8  Integrated  Kernel  Matrix 

Define  the  integrated  kernel  matrix  as 

JCap(t)  =  <7 of  £(rQ,r,)^)^(r')</r, 

and  the  moment  matrix  as 

Map(t)  =  a0  [  M(ra,  r',  t)rl>p(. r')  dr1 

where,  by  definition, 

JC(ra,T',t)  =  J  G(T,r',t)dt 
M(r,r',t)  =  J  g(r,i',t)tdt 

Kkim  CO  =  f  Pum(t)  [-VoOo5pq<f>(r,  t)  +  dpdq${v,  t)]  dr 


(3.47) 


and 

Thus 


and 


m&’m = -  /  x&Sw  * 

In  order  to  evaluate  (3.52)  and  (3.53),  note  that 


47T 


and 


d2x 


=  4 


a3R2 


*(i) -*[*(=)+-(  i)l 


1  - 


3(x  —  a;') 


A2n 


R2 

R\  \  (x  -  x'){y  -  y') 


R 3 


(3.48) 

(3.49) 

(3.50) 

(3.51) 

(3.52) 

(3.53) 

(3.54) 

(3.55) 
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For  the  singular  elements,  we  introduce  a  cubic  exclusion  volume  f2o>  each  side  having 
a  length  2 a  that  is  less  than  or  equal  to  the  smallest  of  the  cell  dimensions,  5X  8y  and 
8Z.  The  exclusion  volume  is  centered  on  the  singularity  at  r  =  0.  Define  the  regularized 
singular  integral  corresponding  to  the  tensor  components  of  Eq.  (3.48)  as[14] 


Kkimit)  =  [  Pkim(r)[-Wo5Pq<t>(T,t)  +  dPdg<!>(r,t)}dr 


(3.56) 


where,  suppressing  the  time  dependence 

Akit  =  [  Pklm(r)dpdq$(r)  dr  (3.57) 


d(P?) 

nklm 

—  j  [/3jdm{r) 

Pklm(0)]dpdMr)dr  +  pklm{ 0)  /  dpdq[$(r)  -  $0(r)]dr 

J  Ho 

(3.58) 

° him 

=  Pklm(Q)  1 9pdq 

[  4>0(r  -  r')  dr' 
Jo  o 

r—O 

(3.59) 

nW 

uhlm 

=  —fJ>QCTQ8pq  I 

Jnkl. 

<£(r )Pklm{r)  dr 

m 

(3.60) 

Here  <J»0(r)  is  the  static  potential  and  $(r)  the  dynamic  potential  defined  as 


$o(r)  = 


47r  |r  | 


*M  = 


erfc(|r|/ o;) 
47r|r| 


(3.61) 


The  contributions  to  Eq.  (3.56)  are  examined  below,  term  by  term  and  case  by  case. 


3.8.1  Evaluation  of 

The  calculation  of  A ^  is  done  using  any  suitable  numerical  quadrature  scheme  which 
approximates  the  integral  over  the  region,  ttkim  —  ^o;  by  a  summation [15].  Thus, 


=  f  n  /W')  W(r)  dr  (3.62) 

»  '52cvPkim(rv)[dPdq®{r)]Tv,  (3.63) 

IV 

where  cv  and  r„  are  the  appropriate  quadrature  weights  and  points  of  evaluation  respec¬ 
tively. 


3.8.2  Evaluation  of 

With  reference  to  Eq.  (3.59),  the  only  case  where  /?Wm( 0)  is  nonzero  is  that  where  k,l,m  = 
0.  Evaluating  the  relevant  integral  and  noting  that  A)oo(O)  =  1.0,  it  is  found  that 


C®  =  A)oo(0)  \dpdq  f  $o(r  ~  r)  dr 


Jr=0 


where  Spq  =  1  if  p  =  q  and  5pq  =  0  otherwise. 


(3.64) 
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3.8.3  Evaluation  of  D^| 

For  a  general  element 


Dkii  =  -Mo<7o<W  /  0M/Wr)  * 


^kln 

=  /o+i,fe  (i  - 1«  _  >i)  (i  - -  *i)  (i  - 1= . -  d) 

J(j-l)Sx  J(k-l)Sy  J(l-l)Sz  \  dx  )  \  dy  /  V  0Z  ) 


U(x)U(y)U(z)  dx  dy  dz 


cr 


^0,(£M§K’(I) 


(3.65) 


where  Fj°\£)  is  given  by  (3.45). 


3.8.4  Evaluation  of 

We  begin  this  section  by  considering  .  Results  for  this  initial  case  are  used  in  the 
evaluation  of  B^  where  k,  l  and  m  take  on  values  that  are  combinations  of  zero  and 
one.  In  all  other  cases  there  is  no  overlap  with  the  exclusion  zone. 

Note  that  the  first  term  in  brackets  in  equation  (3.58)  is 


[/Wr')  -  A)oo(0)] 


s' i  i y1 1  /m x'y' i  ,  i x'z' i  ,  i y'z' i  i x'y'z'  i 

T\ ~  I  a-  I  -  IfI  +  In-'  +  1 77"  I  +  1 77“  I  “  1 777"  I-  (3-66) 

0X  Oy  0z  dXdy  dX0z  dydz 


Substituting  into  the  defining  relationship  for  Bqqq ,  Equation  (3.58),  we  find  that 


?(pp) 

5ooo 


8(-cg?o  -  cfqi  +  c|q2  +  Cgffe  -  c3q4  +  95) 


(3.67) 


where 


cf  =  - 
°  5* 


cx--  +  - 

Cl  ~  5y  +  6, 


dx8y 


+ 


S-nSx 


cl  = 


a* 

W~z 


C3  = 


5X5ydZ 


(3.68) 


and  the  q's  are  defined  by  integrals  given  below.  Note  that  the  integration  in  each  case 
is  over  the  region  defined  by  0  <  x  <  Sx,  0  <  y  <  5y  and  0  <  z  <  5Z,  the  region  being 
referred  to  as  the  first  octant,  fix.  Thus 


1 

qo  =  ~ 

a . 

f  d2$  , 

/  x  x  2  dT 

In  i  ox1 

II 

55  1  l—1 

•  a2#  ,  i  / 

d2$  , 
xy  dr 

ax2 

(3.69) 

O 

oo 

II 

5^^ 

,yzdx2  * 

£ 

II 

d2$  1 

xyzgx,dv  qs  =  l 

r,,  * 

(3.70) 

These  integrals  are  approximated  by  power  series  in  A  =  a/ a,  by  expanding  the  function 
erfc (/?/«)  as  a  power  series  in  R/a  and  integrating  term  by  term.  In  general 


Qs  = 


d2  gjR/ay 

dx2[  R  . 


dr 


(3.71) 
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In  this  form  one  can  compare  similar  developments  for  the  time-  harmonic  case,  where 
1/g!  =  ik  and  g(ikR )  =  exp(ikR),  with  the  transient  case  for  which  g(R/a )  =  erfc (R/a). 
Formally  integrating  with  respect  to  x  gives 


where 


Qa= L  L  y’ z ^ ~  F^0,  V)  z^dy  dz 


F(x,y,z) 


fs( r) 


d  \g{R/a)' 
dx  .  R 


gAM  sWa) 

dx  R 


\df.(x)  ,  *  ,]g(R/a)  ,  x  dg(R/a ) 

vir + &fs{r)\  ~ r~ + MT)w—dir 


(3.72) 


(3.73) 


The  required  series  expansions  are 


g(x)  =  1  -  erf(x)  =  1  -  anx: 

n—0 

where 


Substituting  into  (3.73)  gives 
1 


2n+l 


and 


x 


dgjx ) 
dx 


=  -^an(2n+1)X 

n=0 


(-i)n 


F(x,y,z )  = 


dfjr)  *  1  “r 

Rl  dR  +R?f'^\  +  R  ^ 


VtF  n!(2n  + 1) 

df3(  r) 


n= 0  L 


&  -2n^ 


j?yn+1 


(3.74) 


(3.75) 


(3.76) 


Numerical  integration  of  (3.72)  will  now  give  the  coefficient  in  the  power  series. 

Alternatively  the  integration  can  be  performed  numerically  in  three  dimensions  in 
which  case  we  need  to  evaluate 


cP 
dx 2  L 


g(R/a) 


R 


1  .„  ,  lr,  3x2 

=  -jplRg  -9)1 1 


£ 


_  it 

R?1  R?9 


Using  the  power  series  expansion  for  the  error  function  gives 
dx2 


\g(R/a)' 

1 

’  3x2' 

l  00  r 
- -Y2n 

R  J 

“  ~R? 

.  _  tf2". 

m  ^ 

^  n= 1  L 

Jl  n 


ar  {R\2n+l 


Similarly 


d2 


g{R/a) 


dxdy  [ 

In  this  way  we  get 


R 


3  xy 
R 5 


xy^A  ,  n  (R\2n+1 


(3.77) 


(3.78) 

(3.79) 


q0  «  .035475  +  0.029931A2  —  0.038911A3  +  0.027508A4  —  0.014249A5  (3.80) 

»  —0.0177375  +  .029931A2  —  .032924A3  +  .022662A4  —  0.11850A5  (3.81) 

q2  »  0.0038345  +  .014966A2  -  .020952A3  +  .016141A4  -  0.009015A5  (3.82) 

q3  «  —.007669  +  .014966A2  —  .017959A3  +  .013362A4  —  0.007483A5  (3.83) 

qA  »  .007483A2  -  .011224A3  +  .009353A4  -  0.005612A5  (3.84) 

q5  «  0.059862A2  —  0.059862A3  +  0.037913A4  —  0.018465A5  (3.85) 


23 


]d{xx) 

nooo 

—j^xxT 

-“lOO 

jd{xx) 

■°010 

—j^xxT 

£>001 

■n{xx) 

■°110 

r>{xx) 

■°011 

jd{xx) 

n101 

T){XX) 
n  111 

can 

0 

0 

br> 

8  % 

-4cf 

-4c| 

-44 

2  <4 

0 

2cf 

0 

b(3xx} 

8cf 

0 

-4  cf 

-4  cf 

0 

2c3 

0 

0 

ip- 

— 8c3 

4c3 

4c3 

4c3 

— 2c3 

-2c3 

— 2c3 

C3 

Axx) 

°o 

8 

0 

0 

0 

0 

0 

0 

0 

Table  3.1:  Coefficients  for  the  evaluation  of 


As  the  second  order  partial  with  respect  to  x,  y  or  2  produces  an  even  function  and  /?ooo 
is  even,  contributions  from  the  remaining  7  octants  are  the  same  as  those  for  the  first 
octant.  The  above  considerations  allow  us  to  determine  B^K  In  a  similar  way,  the  other 
possibilities  are  expressed  in  the  form 

Bita  (3-86) 

s=0 

with  the  b's  given  in  Table  3.1.  The  corresponding  results  for  yy  and  zz  components  are 
found  by  appropriate  cyclic  rotation  of  the  coordinate  indices. 

Next  consider  the  evaluation  of  Eq.  (3.58),  for  cases  involving  a  mixed  derivative, 
for  example  derivatives  with  respect  to  x  and  y.  It  is  necessary  to  evaluate  both  of  the 
integrals  in  equation  (3.58).  In  doing  so  one  notes  that  the  mixed  derivatives  are  odd 
functions  in  x  and  y,  therefore  the  first  integral  is  zero  for  all  cases  except  and  . 
The  second  integral  is  zero  for  all  k,  l  and  m,  since  /3fcim(0)  =  0  unless  k  =  l  =  m  =  Oin 
which  case  the  integral  is  zero  by  anti-symmetry.  Hence  we  need  only  concern  ourselves 
with  the  terms 

1  f  d2$  1  f  d2$ 

to  =  -5  /  arj/w-w-  dr  fi  =  /  xyz -5-5-  dr  (3.87) 

a 1  J n,  axoy  a 6  Vn,  oxay 

Expanding  the  integrand  as  a  power  series  in  R/a  and  integrating  term  by  term  gives 


to  »  .020833  -  .007982A3  +  .008741A4  -  0.005714A5  . . . 

ti  «  .006367  —  .003991A3  +  .004846A4  —  0.003448A5 . . .  (3.88) 

Thus  we  express 

(3-89) 

3=0 

with  the  coefficients  are  given  in  Table  3.2. 


3.9  Conclusion 

A  volume  integral  formulation  has  been  presented  for  calculating  the  interaction  of  pulsed  eddy 
currents  with  a  volumetric  flaw  in  a  conductor.  The  approach  uses  a  hypersingular  kernel,  as 
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R  W) 
“100 

r^vJ 

“010 

R(*yJ 

“ooi 

R  W) 
■“110 

R  W) 

“Oil 

R(*y) 

“101 

r>(xy) 

bixy) 

0 

0 

0 

0 

2c| 

0 

0 

0 

b(xy) 

0 

0 

0 

0 

— 2c3 

0 

0 

c3 

Table  3.2:  Coefficients  for  the  evaluation  of 


in  the  corresponding  time-harmonic  problem[16].  This  means  that  the  methods  used  by  Lee  et 
al  [14]  for  the  crucial  calculation  of  the  self-element,  can  be  adopted  for  treating  the  time  domain 
problem[17].  Here  complete  analysis  of  the  self-element  evaluation  has  been  given. 


Appendix  A 
Half- space  matrix 


A1  Half-space  kernel 

In  Sections  7  and  8,  matrix  elements  are  evaluated  from  the  unbounded  domain  integral  kernels. 
In  this  Appendix,  the  matrix  elements  for  half-space  kernels  are  evaluated.  For  a  half-space,  the 
primary  Green’s  function  is  written  as  [9] 

0(r,r  +  &(r1r',tJt')  +  ~Vx  £[V'  x  zV(r,r',t,t')],  (Al) 

where  Go(r,r',t,t')  is  the  free-space  integrated  kernel  given  by  (3.25)  with  (3.11)  and  (3.12). 
The  corresponding  image  term  is  given  by 


Gi(r,r',t,t')  =  - 


/io^'T-VV' 

at  do 


0(r,r",t,t'), 


(A2) 


where  X  —  xx  +  yy  —  and  r"  —  r  —  2 zfz:  is  the  image  point.  In  evaluating  the  Green’s  kernel 
we  note  that  [9] 

V(r,r',t-t')  =  ^y2i  ^V(K,s)e^t~t>'1  ds^j  J0{Kp)KdK  (A3) 

(A4) 


where 


p2  =  (x  -  x ')2  +  {y  -  y1)2  and  7  =  —  po<tqs. 


-l\z+z’\ 


A2  Moment  kernels 

The  integrated  kernel  is  defined  as 

/C(r,  r',  t,  t?)  =  J  Q(r,r 1  d(t  - 11) 

and  the  (first  order)  moment  kernel  as 

M(r,r',t,t')  —  J(t-  t')G(r,r',t,t')d(t  -  t'). 
Therefore  the  integrated  kernel  is  given  by 


(A5) 


(A6) 


fC(r,  r',  t,  t')  =  fCo(r,  r',  t,  t')  +  /Cj(r,  r',  t,  t')  +  — V  x  £[V'  x  zV (r,  r',  t,  t%  (A7) 

o’o 

where  K,o(r, is  the  free-space  integrated  kernel  given  by  (3.24)  with  (3.11)  and  (3.12). 
The  corresponding  image  term  is  given  by 


/Cj(r,  r',t,  t')  =  - 


p0<f>(r,  r ",  t ,  t')X'  +  — W'$(r,  r",  t,  t') 
<?Q 


(A8) 
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The  first  moment  kernel  is  given  by 


M{r,r',t,t')  =  Mo{r,r'  +  Mi(r,r',t,t')  + 


—Vxz[V'xzW(r,r',t,t% 

&0 


(A9) 


where 


dt 

roc 


(27 r)2z  Jo 
(WO  f°°  f  _1_ 
(27r)2i  Jo  Jbt  27 s 


1  dV(n,  s ) 


9s 


7 


s*(*  ^  ds  j  Jo(kp)k9k 

<M) 


ea(t-t')_7|z+z'|  ds  JQ(Kp)KdK 

(AlO) 


A3  Integrated  Kernel  Matrix  evaluation 

The  matrix  for  the  integrated  kernel  is  written  as 


ICklm,KLM(t )  —  +  ^k-K,l-L,m+M+  lCO  +  ^-KV-L.m+M+lW) 

where  A:,  K  =  0,  1,  2 . . .  7VS,  /,  L  =  0,  1,  2 . . .  Ny,  m,  M  =  0,  1,  2,...NZ.  4-Ar,i-L,m-M(*)  is 
given  by  (3.52)  and 


-(0 


(V) 


K$H(t)  =  --  f  Pklm( r') V'  x  z[V  x  zV( 0,  r',  t, 0)]  dr' 

=  —  r+l  r  r+1  dx'dy'dz'pjf )(h(f \pm[y )v'  x  z[V'  x  zV(0, r', t, 0)] 
(Tq  Jxk-i  Jyi~\  Jzm- 1  \0x/  /  \^/ 


f07  [ 

yx}klm V 


-KS3u»(‘)  0 


0 


0  0 


(A12) 


and  we  have  written  =  m<5z  —  2d,  zm- 1  =  zm  —  Sz  and  zm+i  =  4-  5Z. 

Thus,  the  matrix  contains  only  transverse  components  and  can  be  written  as 


=  ^SIZ  C  C  dzdxdy^{i)  *  (t)  ^(£) 


z\  d2V (0, r, t, 0) 


dpdq 


where  p,  q  =  x  or  y. 
Note  that 


72  K2n2  k 

Of  iJoM}  =  -  - 


2p2 


Ji{np) 


,  (A13) 


(A14) 
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and 


which  suggests  that  we  can  write 


Krhmit)  -  ao  JXki  Jyii  if  J  S°  (P> m ’  *)  +  f 1  p2  J  ^ m’  *) 


j9fc(?)  Ajf-Uyda; 
\()/  \°y  J 

(A16) 


=  ~J  ^  f‘+1  (p)  [5S(p,m,t)  -25}(p,m,t) 

u0  ''xjt-i  •'J//-X  \p  / 


£ )  ft (|-  ]  ds  dy.  (A17) 

v.  /  \  y 


“  <2^i  r  l  j g  -  /m(,s' s)e*' <,a  <ai8) 


where 


rZm+ 1  /  z  \ 

/m(«,a)  =  /  A»  I  jS-  exp(-7*)<fe 

J  Zm-\  \°Z/ 


~2  exP(— 7zm)[2  -  exp  (7^)  -  exp(— 7<iz)] 


(A19) 


Write 


«1/^’  5  =  [  [2*8  (a, «,  *m)  -  ^3 (a,  «,  *m+l)  -  ^3 (a,  K,  - 


kJq(kp) 

l/pJx(«p) 


where  (see  Appendix  B) 


(A20) 


(A21) 
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A4  Moment  Kernel  Matrix  evaluation 


pklm(r')V  x  z[V  x  zW(0,r',tMdr' 

a0  J^klm 

=  —  r+i  r  f m+1  dx>  dyi  dzt^  {  t )  &  (  f W  (  f) v' x  ^v' x  * w{o’ r>’ ^ °)] 

&0  Jxk-i  Jyi- 1  Jzm- 1  \dx/  \ /  V^/ 

-*£*»<*)  "Skw  ° 


-J^jarnW  0 

0  0  0 

Thus,  the  matrix  contains  only  transverse  components  and  can  be  written  as 

C!  C izdxdyh  (I) (t) **  (i)  (A23) 


(A22) 


where p:q  =  x  ov  y.  We  write 
1  r®fc+i  rw+i 

"£&.<*)  =  --/  / 

°0  Jxk- 1  ^y;_i 


^  -?o(p,m,t)  +  f  1  -  5i(p,m,i) 


■(f) 


A I  7  )  A I  ^  J  dv  dx 

(A24) 


where 


-  —  f  k+1  f m+l  [fig(p,m,t)  —  2S?(p,m,t)]  da:dy-  (A25) 

0“O  Jxk~i  Jyi- 1  \P  /  ^  J  \ ^ac /  \®»/ 


502(p,m,t)  1  =_J_  r  f  f  (K  Aestds 

sf(p,m,t)  (27r)2i  70  ABr  \  s/5s\Vk  7/  ’  J 


k2J0(kp) 

K/pJi(np) 


ndn  (A26) 


Write 


S§(p,m,t ) 
£%{p,m,t) 


=  2tt6~  Jo  \.2F(a’K’z™)  ~F(a,K>zm+l) -F(ot,K,Zm-i)] 


where 
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Appendix  B 
Half-space  matrix 


B1  Definitions  of  Functions 

For  convenience,  have  defined  a  function  i^(a, £)  through  its  inverse  Laplace  transform.  Thus 

(B1) 

which  can  be  written  as 

Fv(a,  k,  0  =  e-«2(*-t')//w o  J_  f  L - e-V*<+»(t-0/«w>dg 

27TZ  JBr  S1'/2 ( +  «) 

For  example,  it  is  found  that  [15] 


Fi(a,n,0  =eKCerfcf^  +  ^V 


B2  Recursive  Relation 

Note  that 


1  1 


7^(7  +  «)  K 


1  1 


Y  7*'"1(7  +  «)J 


therefore 


Hence 


Fv(a,K, C)  = 


1 


exp  (— a2K2/4)a"  2iv  2erfc(^)  —  Fl/_i(a,  k, 


F2(a,K,C)  =  - 

K 


which  can  also  be  written  as 


F2(a,K,  C)  =  — e  52 

K 


=.C2/«2erfc  ^  —  e^2erfc(£)  , 


(B2) 


(B3) 


(B4) 


(B5) 


-a2K2/4erfc^^  _  cleric  ^  +  “)  ’  (B6) 


(B7) 


/•  2 

where  £  =  *  +  Qr-  Note  that  ex  erfc(z)  may  be  evaluated  at  large  arguments  using  an  asymptotic 
expansion[15]. 


B3  Further  Relationships 


Define 


and  note  that 


««■*■<> =hL$-rM,-e)is 


1  fiocr  I"  1  1 

s  2k  Lt  —  «  7  4-  k\ 


(B8) 

(B9) 
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hence 


fu(a,K,0 


f  J_[_J _ L_ 

2 k  2m  Jbt  'f  1 7  —  «  7  +  k. 

-K,()  «,0] 


e-rt+s(t-t‘)ds 
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Veetrans  Installation  Guide 

Veetrans  is  a  computer  code  for  calculating  the  pulsed  eddy  current  signals  due  to  a  region  of  corrosion 
in  an  otherwise  uniform  conductor.  The  signals  are  measured  as  a  function  of  time  by  a  Hall  sensor  on 
the  axis  of  an  excitation  coil. 


Code  Installation 


1 .  Initially  the  files  should  be  copied  into  one  folder  CWeetrans. 

2.  Uses  Windows  Explorer  or  My  Computer  to  find  install.bat  in  C:\Vtrans  and  double  click  on  the 
icon  with  the  mouse. 

3.  The  install.bat  will  rearrange  the  files  putting  the  main  executable  VtransGMats.exe  in 
C:\Veetrans\bin  along  with  other  files  associated  with  the  execution. 

The  application  can  be  started  via  Windows  Explorer  or  My  Computer  or  by  adding  VtransGMats  to  the 
Task  Bar  or  the  Desk  Top  in  the  usual  way. 


Folders  Structure 


1 .  By  default,  the  path  to  the  root  folder  is  C:\Veetrans. 

2.  The  folder,  Veetrans,  has  two  sub-folders,  bin  and  tips. 

3.  The  folder,  bin,  contains  the  application  which  consists  of  an  executable  file,  VtransGMats.exe  and 
other  files  supporting  the  execution.  In  the  folder  tips  files  containing  input  and  output  data  are 
stored.  The  main  input  files  define  the  probe,  the  probe  scan  and  the  flaw.  These  files  have  the 
extension  .tip  indicating  that  it  is  a  transient  data  input  file.  The  main  output  file  has  the  extension 
.tm  indication  that  is  contains  the  coordinates  of  the  probe  and  transient  signal  data. 

A  .tip  is  edited  to  get  the  desired  input  data  and  the  output  data  found  in  the  .tm  file  can  be  plotted  using 

suitable  graphics  utilities  (not  supplied). 
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Veetrans  User  Guide 

This  guide  describes  how  to  use  Veetrans  to  compute  the  transient  eddy  current  probe  signal  due  to  a 
region  of  corrosion.  The  signals  are  measured  as  a  function  of  time  by  a  Hall  sensor  on  the  axis  of  an 
excitation  coil.  The  output  data  are  magnetic  field  values  in  Amps/m. 


Veetrans  Calculations 


The  calculation  of  the  transient  probe  signal  due  to  a  region  of  corrosion  depends  on  (a)  the  probe 
parameters  (b)  the  specimen  parameters  (c)  the  position  of  the  probe.  These  parameters  are  specified 
in  a  .tip  file.  Several  such  files  can  be  created  and  the  calculations  for  each  file,  called  a  job,  carried  out 
in  sequence.  A  sequence  of  jobs  is  called  a  run.  The  .tip  files  used  in  for  a  run  are  specified  in  the 
Runspec.ctr  file  located  in  the  tip  folder.  This  is  called  the  control  file  because  it  controls  the  selection  of 
the  .tip  files  that  are  used  in  a  particular  run  (sequence  of  jobs).  Below  is  an  example  of  a  Runspec.ctr 
file. 


# 

Number  of  jobs  (integer) 

2 

Job  number  (integer)  File  root  name 

0  'circlOxlOxl' 

1  'Slotlxl0x6' 


« Figure  1 :  A  typical  Runspec.ctrfWe 

The  above  control  file  will  take  the  data  from  the  file  named  circlOxlOxl. tip  and  performs  the 
calculation  defined  in  it.  It  then  does  the  same  thing  for  the  file  Slot1x10x6.tip.  Note  that  you  are 
required  to  specify  the  number  of  jobs  to  be  performed,  a  job  number  for  each  and  the  file  root  name  for 
the  job.  Notes  can  be  added  to  the  file  prior  to  the  #  and  these  will  be  ignored  by  the  application.  This 
allow  the  user  to  write  some  comments  in  the  file.  There  must  be  a  file  circlOxlOxl. tip  and  another 
Slot1x10x6. tip  in  the  tips  folder  otherwise  the  application  will  not  be  able  to  perform  the  calculation.  It  is 
also  important  that  the  .tip  files  have  the  correct  format. 
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Before  a  calculation  is  performed  it  is  usually  necessary  to  create  or  edit  a  .tip  file  to  specify  the 
calculation.  An  example  of  such  a  file  is  given  below.  The  application  begins  to  acquire  data  from  the 
file  after  the  #  0  is  encountered.  Before  that  point  any  information  can  be  written  into  the  file  but  it  is  not 
used  by  Veetrans.  There  follows  eight  data  fields  labelled  0  though  7. 


#0 

ONOTE 

1x10x6  flaw  grid,  flaw:  0.25x9. 97x8. 016mm.  30  x  20  microsec. 

1  TIME 

Sat  Jun  24  2000 

2  PROBLEM  Root  problem  domain  accuracy  conductivity(S/m) 

'slot  10x6'  1  1  0.01  1.8792E+07 

3  FLAW  no  x  no  y  no  z  dx  (mm)  dy  (mm)  dz  (mm)  setflaw 

1  10  6  0.25  0.997  1.336  -1 

100  100  100  100  100  100  100  100  100  100 

100  100  100  100  100  100  100  100  100  100 

100  100  100  100  100  100  100  100  100  100 

100  100  100  100  100  100  100  100  100  100 

100  100  100  100  100  100  100  100  100  100 

100  100  100  100  100  100  100  100  100  100 

4  SCAN  startx(mm)  starty  (mm)  step  x  (mm)  step  y  (mm)  no.x  no.y  lift  off  (mm) 

0.0  0.0  2.5  2.5  5  5  0.6 

5  PROBE  COIL  inner  r(mm)  outer  r  bottom  z  top  z  no.  turns 

7.1  11.41  0.0  5.0  2550 


0  NOTE:  This  is  a  one-line  text  string  acting  as  a  summary  of  the  file  data. 
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1  TIME:  The  date  and  time  of  the  file  creation 


2  PROBLEM:  This  data  field  contains: 

1 .  The  root  name  of  the  output  files. 

2.  An  integer  designation  of  the  problem  type.  Always  use  1  to  indicate  that  a  forward  problem  is 
being  solved. 

3.  The  domain  number  is  also  always  1.  This  denotes  that  the  conductor  is  modelled  as  a  half¬ 
space. 

4.  The  accuracy  is  then  given  as  a  decimal  fraction.  A  value  of  0.001  is  suitable. 

5.  The  conductivity  of  the  material  is  then  given  (in  S/m). 

3  FLAW:  The  flaw  is  made  up  of  a  number  of  cells  in  a  regular  array.  Each  is  a  rectangular 
parallelpiped  of  the  same  size.  The  flaw  data  field  defines  the  number  of  cells  in  the  x-,  y-,  and  z- 
directions  and  the  size  of  the  x-,  y-  and  z-dimension  in  mm.  The  value  of  setflaw  is  always  -1.  The  z 
coordinate  is  the  vertical  direction. 

The  conductivity  of  each  ceil  is  denoted  by  an  integer  between  1  and  100.  A  value  of  1  means  that  the 
cell  has  a  conductivity  1  %  less  than  that  of  the  host  material,  a  value  of  100  means  that  the  conductivity 
of  the  cell  is  100%  less  than  that  of  the  host  and  therefore  zero,  a  value  of  75  means  that  the  cell 
conductivity  is  25%  of  the  host  value. 

In  the  example  above  all  the  cells  have  a  100%  reduction  in  conductivity  and  hence  they  have  zero 
conductivity.  This  is  typical  of  the  way  a  crack  or  a  void  due  to  corrosion  is  represented  since  they  are 
usually  assumed  to  have  zero  conductivity.  The  cell  values  are  given  in  blocks  separated  by  a  space. 
Each  block  has  the  same  number  of  cells  and  represents  different  levels  in  the  flaw  starting  at  the  top 
and  working  down.  In  the  example  there  is  only  one  line  of  values  in  each  block  indicating  that  there  is 
only  one  cell  in  the  x-direction.  There  are  10  values  in  the  line  one  for  each  of  the  cells  in  the  y-direction 
and  there  are  6  data  blocks  indicating  that  there  are  6  levels  of  cells  in  the  z-direction.  The  example 
therefore  represents  a  slot  of  width  0.25  mm  in  the  x-direction,  9.97  mm  long  and  8.016  mm  deep. 

Please  note  that  none  of  the  cell  values  can  be  zero.  This  will  cause  problems  with  the  calculation. 

4  SCAN:  The  probe  position  is  varied  over  a  regular  grid  starting  at  a  location  defined  by  startx  and 
starty.  These  are  the  coordinates  of  the  probe  center  in  the  xy-plane.  The  step  size  in  the  x  and  y 
directions  are  given  in  mm  followed  by  the  number  of  steps  in  the  x  and  y  directions.  Finally  the  probe 
lift-off  is  defined.  This  is  the  distance  from  the  coil  reference  level  to  the  top  of  the  work  piece. 

5  PROBE:  The  probe  consists  of  a  coil  with  a  hall  sensor  on  its  axis.  In  the  probe  data  field,  the  inner 
and  outer  radius  of  the  coil  is  defined,  the  z  coordinate  of  the  bottom  and  the  top  of  the  coil  and  the 
number  of  turns.  Note  that  the  lift-off  value  is  added  to  the  top  and  bottom  coordinates  to  get  the 
location  of  the  top  and  bottom  of  the  coil  with  respect  to  the  surface  of  the  work  piece. 

6  TIME  STEPS:  This  defines  the  number  of  time  steps  used  to  describe  the  pulse  and  the  time  interval 
between  each  signal  sample.  This  time  interval  must  be  carefully  controlled  as  discussed  in  the  section 
Getting  Started. 

7  FLAW  DEPTH:  The  first  value  in  this  data  field  is  the  distance  from  the  top  of  the  flaw  region  to  the 
surface  of  the  conductor.  This  distance  is  designated  the  ligament.  The  second  value  is  the  rise  time  of 
the  coil  current.  The  excitation  of  the  pulsed  eddy  current  is  due  the  changing  field  caused  by  the 
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variation  of  the  current  in  the  coil.  It  is  assumed  that  the  coil  current  increases  exponentially  in  time  with 
a  time  constant  given  in  this  data  field.  Typically  values  for  aircraft  inspection  for  corrosion  or 
subsurface  cracks  are  in  the  range  50-500  micro-seconds.  The  third  value  in  this  data  field  is  the  height 
of  the  Hall  sensor  above  the  work  piece  in  mm  when  the  liftoff  is  zero.  The  Liftoff  gets  added  to  the 
height  of  the  Hall  sensor  to  give  the  distance  from  sensor  to  work  piece. 


Getting  Started 


The  basic  procedure  for  running  a  calculations  is  as  follows: 

1 .  Create  or  preferably  edit  a . tip  file  to  define  the  input  data. 

2.  Edit  the  Runspec.ctr  file  to  include  the  new  data  file  in  the  list  of  jobs. 

3.  Run  the  application. 

The  calculation  becomes  more  accurate  if  the  cell  size  is  small  but  the  smaller  the  cell  size,  the  greater 
the  number  of  cells  needed  for  a  given  flaw  and  the  longer  it  takes  to  do  the  calculation.  It  is  best  to  do 
small  calculation  until  experience  with  the  results  has  been  gained.  Large  calculations  may  have  to  run 
overnight  but  it  is  best  to  have  some  experience  running  small  jobs  before  making  the  commitment  to 
large  calculations. 

The  calculation  becomes  more  accurate  if  a  small  time  interval  is  chosen  but  again  the  increase  in 
accuracy  comes  at  the  cost  of  a  longer  calculation.  A  good  rule  of  thumb  is  to  assign  a  time  increment 
dependent  on  the  largest  cell  dimension:  let  da  =  max(dx,  dy,  dz)  then  the  time  increment  is 
approximated  by  dt=p.a  (da)2. 
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Consulting  Set  vices 

I.  INTRODUCTION: 

The  capability  of  a  nondestructive  inspection  is  dependent  on 

•  Flaw  (Artifact)  Variables 

•  Test  Object  Variables 

•  NDE  Method  Variables 

•  NDE  Materials  Variables 

•  NDE  Equipment  Variables 

•  NDE  Procedure  Variables 

•  NDE  Process  Variables 

•  Calibration  Variables 

•  Acceptance  Criteria  /  Decision  Variables 

•  Human  Factors 

A  nondestructive  inspection  measurement  is  not  single  valued,  but  is  probabilistic  in  nature  due  to  its 
multi  -  parameter  nature.  The  basis  for  use  of  probability  of  detection  (POD)  as  a  metric  is  discussed 
in  detail  in  the  attached  papers  and  will  not  be  discussed  here  [REF  1,2]. 

Since  POD  was  developed  and  has  been  extensively  applied  to  the  assessment  of  crack  detection 
capabilities  in  advanced  engineering  structures,  examples  of  use  in  crack  detection  are  included  in  the 
example  data  sets. 

NOTE:  Copies  of  the  references  are  included  on  the  CD. 

II.  REQUIREMENTS  AND  EXAMPLES  INCLUDED: 

A.  SETTING  UP  THE  COMPUTER  FOR  DATA  ANALYSIS: 

Excel  ™  95/7.0;  97  and  2000  have  built  in  charting  capabilities  and  templates.  It  also  provides  for 
user  templates  when  a  custom  analysis  and  plotting  routine  is  desired.  Chart  templates  must  be  loaded 
BEFORE  the  analysis  process  is  started.  Set  up  instructions  and  templates  for  both  shell  and  crack 
analysis  are  included  on  the  CD  (File  -  Computer  Set  up. doc).  The  following  page  is  a  copy  of  the 
single  sheet,  set  up  procedure. 
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SETTING  UP  A  COMPUTER  FOR  ANALYSIS 

Excel  programs  have  built  in  charting  capabilities  and  templates.  They  also  provides  for  user 
templates  when  a  custom  analysis  and  plotting  routine  is  desired.  User  templates  must  be  loaded 
BEFORE  the  analysis  process  is  started. 


To  load  a  user  template: 

a)  Select  an  analyzed  data  file  with  a  plot  in  the  format  desired.  Template  plots  have  been  provided 
for  analysis  of  crack  data  by  crack  length.  Load  Excel  and  the  appropriate  chart  example.  The 
chart  in  the  proper  format  must  be  visible  when  setting  up  a  template. 

b)  From  Tool  Bar  >  [FORMAT]  or  [CHART] 

c)  Select  [AutoFormat]  or  [Chart  Type] 

d)  Select  [User  Defined]— 

e)  Select  [Customize] 

f)  Select  [Add]  - 

g)  Dialog  box  will  appear  -  INPUT  [chart  name  —  for  example  -  UTJLength 
Use  this  user  format  in  subsequent  plotting. 
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B.  FILES  INCLUDED  ON  THE  CD 

1 .  ®  [Sample  Plots] 

Penetrant  lA.xlm 
Penetrant  lB.xlm 
Penetrant  lC.xlm 
Penetrant  2A.xlm 
Penetrant  2B.xlm 
Penetrant  2C.xlm 
Penetrant  3A.xlm 
Penetrant  3B.xlm 
Penetrant  3C.xlm 


2.  ■  COMPUTER  SET  UP 

W  -  Computer  Set  Up.doc 

3.  □  Manual  in  text  format 

W  -  POD  Instructions 

W  -  Reference  :  Rummel  -  Plenary  Honor  Lecture  Paper  -  1999.doc 

4.  ■  Crack  Length  Analysis  (The  MACRO  TOOLS) 

XL  -  Interpol.xls  -  An  interpolation  tool 
XL  -  Maxlikel.xls.  -  Maximum  likelihood  analysis  module 
XL  -  Verify.xls  -  Data  verification  module 

5.  ®  Sample  Files 

XL  -  Sample  A  -  Order.xls 
XL  -  Sample  B  -  Order.xls 

6. °  [POD  INSTRUCTIONS] 

Excel  95.7.0 
Excel  97 
Excel  2000 
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e.  At  the  top  of  the  columns,  Select  Columns  A,  B,  C,  D,  and  E.  (They  will  appear 
highlighted) 

f.  On  the  TOOL  BAR  select  Z-A-l.  All  data  will  now  be  ordered  in  concert  to  the  data  in 
Column  A. 

g.  Click  on  Column  A;  right  click  and  select  DELETE.  Data  will  now  be  in  the  form 
required  for  Analysis. 

For  your  convenience  and  exercise,  unordered  data  is  provided  in  the  file  folder  named 
[Example  Data  Files].  Select  one  of  the  data  sets  in  this  folder  for  practice  or  you  may  prepare 
your  own  files  from  reference  data  or  user  generated  data. 
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2.  REQUIRED  FILES 

For  execution,  the  macro  requires  that  all  three  files  provided  in  the  [Crack  Length 
Analysis]  folders  be  loaded.  It  is  necessary  to  close  these  files  (do  not  save)  each  time  the  Macro 
program  is  run. 

Select  the  folder  [Crack  Length  Analysis];  the  folder  will  open  to  show  the  three  required  files  for 
analysis.  Hold  down  the  shift  key  and  select  all  three  files.  The  files  will  now  be  highlighted.  Select 
OPEN.  The  three  files  will  now  be  loaded. 

Before  Macro  execution,  it  is  necessary  to  have  the  data  file  in  view.  From  the  <TOOL  BAR>,  select 
<WINDOW>  and  click  on  the  name  of  the  prepared  data  file.  Position  the  cursor  at  the  top  of  Column 
E. 

“Ctrl  M”  Starts  the  analysis. 

Follow  the  step-by-step  instructions  for  plotting  as  shown  on  the  POD  INSTRUCTION  sheets. 

The  first  notice  that  may  appear,  is  a  message  that  informs  the  user  that  the  data  distribution  is  not  in  a 
form  that  will  allow  the  maximum  likelihood  analysis  to  converge.  This  indicates  that  the  flaw 
distribution  does  not  meet  the  requirements  for  absolute  operation  of  the  model.  The  data  for  the 
larger  flaws  will  then  be  used  for  the  analysis.  This  will  result  in  a  conservative  POD  plot.  There  are 
several  reasons  why  data  distribution  will  not  meet  the  model  requirements,  but  the  most  common  are 
data  out  of  order  and  too  few  misses  at  the  lower  defect  sizes.  The  distribution  must  be  adjusted  for 
absolute  convergence,  but  may  not  be  possible  in  all  cases.  It  is  advised  that  a  note  be  added  to  the 
footer  information  on  all  plots  that  do  not  converge.  This  lets  the  user  know  that  the  resulting  plot  is 
conservative  and  the  actual  capability  is  better  than  indicated  by  the  data. 

The  second  notice  that  the  operator  must  deal  with  is  the  plotting  parameters  (noted  as  crack  length). 
The  crack  length  analysis  program  has  been  pre-set  to  plot  from  0  to  0.750  inch  in  increments  of 
0.005  inch.  Click  ENTER  to  accept  the  pre-set  values  or  modify  the  values  to  the  desired  plot.  It  is 
often  desirable  to  plot  as  a  function  of  crack  depth  or  crack  area.  The  range  can  be  adjusted  to  the 
desired  values.  It  may  be  convenient  to  save  the  files  in  a  separate  folder  for  crack  depth  or  area 
analyses.  The  analysis  programs  have  been  protected  with  a  password  to  avoid  contamination  and 
overwriting  for  different  applications. 

The  analyzed  parameters  spread  sheet  will  appear  and  the  plotting  routine  should  be  followed  in 
accordance  with  the  POD  INSTRUCTION  sheets. 
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Program  Instructions 
(Applicable  to  Excel  95/7.0) 

1 .  [USER  FILE  NAME]  -Data  must  be  in  the  form  -  Crack  No.  -  Crack  Length  -  Crack  Depth  - 1  or 
0  (Detect  /  nondetect)  -  for  analysis  by  crack  length. 

2.  From  the  <FILE>//<OPEN>  Select  [Crack  Length  Analysis]  folder  and  load  all  three  programs  - 
“verify.xlm”;  “maxlikel.xlm;  and  ‘interpol.xlm.  (for  length).»>{{OPEN}}.  A  Dialog  box  will 
appear.  Click  [READ  ONLY]. 

3.  Bring  data  to  front  (From  the  <TOOL  BAR>  Select  <WINDOW>,  [USER  FILE  NAME],  Place 
cursor  at  the  top  of  the  first  open  column  -  right  of  “Verified”  (Column  E). 

4.  CTRL  -  M;  to  start  analysis.  Analysis  is  completed. 

5.  A  Dialog  Box  appears  “POD  Curve  Table”.  Select  the  “maximum  crack  length  “  and  “crack 
increment”  to  be  used  in  plotting.  0.750  and  0.005  are  the  defaults. 

6.  From  the  <TOOL  BAR>  Click  chart  wizard;  highlight  columns  $L  and  $M 

7.  Click  Insert  -  Chart  -  As  a  new  page. 

8.  A  Dialog  box  appears,  with  range  $L:$M,  Click  next. 

9.  Chart  type  appears  -  click  scatter 

10.  Chart  type  appears  -  click  line  fit. 

11.  Click  finish  -  POD  curve  appears. 

12.  From  the  <TOOL  BAR>  Click  -  insert  new  data  -  a  box  appears  with  range. 

13.  Switch  to  sheet  1,  [Bottom  tool  bar]  select  $0:$P,  click  chart. 

14.  Click,  x  “Categories  (X  values)  in  first  column).Then  OK 

15.  From  the  <TOOL  BAR>  Select  format  -  Autoformat. 

16.  Select,  user  defined 

17.  Click  previous  [USER  FILE  PLOT}  POD  -  OK.  (Must  have  loaded  the  format  previously  -  To 
load  the  applicable  mode-  See  Note  1 . 

18.  Change  descriptive  box  data  as  applicable. 

19.  Enter  90%  POD  value  -  See  Note  2. 

20.  From  the  <TOOL  BAR>  Select  <FILE>  Click  Page  set-up,  Footer,  enter  description. 

21.  Data  analysis  is  complete. 

22.  Print  as  desired. 

Note  1: 

To  load  a  chart  template  (For  Example  -  Pen  1  A): 

a)  Load  a  previous  POD  plot  into  Excel. 

b)  From  Tool  Bar  >  Select  [FORMAT] 

c)  Select  [AutoFormat] 

d)  Select  [User  Defined]—  upper  left  comer 

e)  Select  [Customize]  -  right  comer 
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f)  Select  [Add]  -  right  side 

g)  Dialog  box  will  appear  -  INPUT  [chart  name  —  [FORMAT]  for  example  -  UT_Length 

h)  Select  -  [OK],  Dialog  Box  will  close 

i)  Use  this  user  format  in  subsequent  plotting. 

Note  2: 

To  Interpolate  for  the  90%  POD  Value 

a. )  Got  to  Sheet  1  [  Data  ]  --  on  the  bottom  tool  bar. 

b. )  Scroll  down  Column  M  to  the  values  at  89. XXX;  90.XXX 

c. )  Select  these  four  cells  <Example> 


s 
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M 

0.39 

89.86499 

0.395 

90.01687 

d. )  Copy  these  four  cells. 

e. )  From  the  <TOOL  BAR>  Select  <WINDOW>,  Select  interpol.xls 

f. )  From  the  <TOOL  BAR>  Select  <EDIT>  “Paste  Special” 

g. )  A  Dialog  box  appears.  Select  “Values” 

h. )  The  interpolation  is  completed  and  the  90%  POD  value  appears  to  the  right  in  units  of  both  inches 

and  millimeters. 

i. )  Copy  these  values  to  the  descriptive  data  box  on  the  <Chart>. 

Your  final  Chart  should  look  like  this: 


ACTUAL  CRACK  LENGTH  -  (Inch) 
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Program  Instructions 
(Applicable  to  Excel  97 ) 

1 .  [USER  FILE  NAME]  -Data  must  be  in  the  form  -  Crack  No.  -  Crack  Length  -  Crack 
Depth  - 1  or  0  (Detect  /  nondetect)  -  for  analysis  by  crack  length.. 

2.  From  the  <FILE>//<OPEN>  Select  [Crack  Length  Analysis]  folder  and  load  all  three 
programs  — ‘ “verify.xlm”;  “maxlikel.xlm;  and  ‘interpol.xlm.  (for 
length).»>{{OPEN} }.  A  Dialog  box  will  appear.  Click  [Enable  Macro’s}  —Click 
[READONLY], 

3.  Bring  data  to  front  (From  the  <TOOL  BAR>  Select  <WINDOW>,  [USER  FILE 
NAME],  Place  cursor  at  the  top  of  the  first  open  column  -  right  of  “Verified” 
(Column  E). 

4.  CTRL  -  M;  to  start  analysis.  Analysis  is  completed. 

5.  A  Dialog  Box  appears  “POD  Curve  Table”.  Select  the  “maximum  crack  length  “  and 
“crack  increment”  to  be  used  in  plotting.  0.750  and  0.005  are  the  defaults. 

6.  Highlight  columns  $L  and  $M 

7.  From  the  <TOOL  BAR>  Select  [Insert  Chart} 

8.  A  Dialog  Box  Appears  -  Select  [xy  scatter];  NEXT 

9.  From  the  top  bar  Select  [Series]  NEXT . NEXT 

10.  A  dialog  box  appears.  -  Select  “Place  Chart”  [As  a  new  sheet]. 

11.  Click  [Finish]-  POD  curve  appears. 

12.  From  the  <TOOL  BAR>  Select  Chart  -  Add  Data 

13.  Click  -  insert  new  data  -  a  dialog  box  appears  with  range. 

14.  Switch  to  sheet  1,  [Bottom  tool  bar]  select  SO:$P  -  [OK] 

15.  A  dialog  box  appears  {PASTE  SPECIAL] 

16.  Click,  x  “Categories  (X  values)  in  first  column). Then  OK 

17.  From  the  <TOOL  BAR>  Select  <CHART>  , CHART  TYPE> 

18.  A  dialog  box  appears.  Select,  “user  defined” 

19.  Click  previous  [USER  FILE  PLOT}  POD  -  OK.  (Must  have  loaded  the  format 
previously  -  To  load  the  applicable  mode-  See  Note  1 . 

20.  X’s  appear  on  the  top  and  bottom  of  the  chart  (Hit  /  Miss  Values).  Select  one  of  the 
X’s.;  Double  Click. 

21.  A  dialog  box  appears.  Select: 

22.  Line  should  be  NONE 

23.  Marker  should  be  CUSTOM 

24.  Select  Style  Symbol  (ie.  A) 

25.  Foreground  -  Select  BLACK 

26.  Background  -  Select  BLACK 

27.  Size  -  Select  8  pt. 

28.  [OK] 

29.  If  you  wish  to  have  an  insert  descriptive  box,  Select  from  another  chart  (be  sure  the 
outline  on  the  box  is  fully  cross  hatched.  Copy  and  paste  to  the  new  chart.  Text  can 
then  be  edited  to  annotate  the  new  chart. 
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Program  Instructions 
(Applicable  to  Excel  97) 

30.  Enter  90%  POD  value  -  See  Note  2. 

3 1 .  From  the  <TOOL  BAR>  Select  <FILE>  Click  Page  set-up,  Footer,  enter  description. 

32.  Data  analysis  is  complete. 

33.  Print  as  desired. 

Note  1: 

To  load  a  chart  template  (For  Example  -  Pen  1  A): 

a)  Load  a  previous  POD  plot  into  Excel. 

b)  From  Tool  Bar  >  Select  [CF1ART] 

c)  Select  [Chart  Type] 

d)  Select  [Custom  Type]  -  Top  Bar 

e)  Select  [User  Defined]--  Lower 

f)  Select  [Add] 

g)  Dialog  box  will  appear  -  INPUT  [chart  name  —  for  example  -  UT_Length 

h)  Select  -  [OK],  Dialog  Box  will  close 

i)  Use  this  user  format  in  subsequent  plotting. 
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Note  2: 

a)  To  Interpolate  for  the  90%  POD  Value 

b)  Got  to  Sheet  1  [  Data  ]  —  on  the  bottom  tool  bar. 

c)  Scroll  down  Column  M  to  the  values  at  89.XXX;  90.XXX 

d)  Select  these  four  cells  <Example> 


L 

M 

0.39 

89.86499 

0.395 

90.01687 

e)  Copy  these  four  cells. 

I)  From  the  <TOOL  BAR>  Select  <WINDOW>,  Select  interpol.xls 

g)  From  the  <TOOL  BAR>  Select  <EDIT>  “Paste  Special” 

h)  A  Dialog  box  appears.  Select  “Values” 

i)  The  interpolation  is  completed  and  the  90%  POD  value  appears  to  the  right  in  units  of 
both  inches  and  millimeters. 

j)  Copy  these  values  to  the  descriptive  data  box  on  the  <Chart>. 
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Program  Instructions 
(Applicable  to  Excel  2000) 

1 .  [USER  FILE  NAME]  -Data  must  be  in  the  form  -  Crack  No.  -  Crack  Length  -  Crack 
Depth  - 1  or  0  (Detect  /  nondetect)  -  for  analysis  by  crack  length.. 

2.  From  the  <FILE>//<OPEN>  Select  [Crack  Length  Analysis]  folder  and  load  all  three 
programs  — “verify.xlm”;  “maxlikel.xlm;  and  ‘interpol.xlm.  (for 

length). »>{ {OPEN}}.  A  Dialog  box  will  appear.  Click  [Enable  Macro’s}  —Click 
[READ  ONLY], 

3.  Bring  data  to  front  (From  the  <TOOL  BAR>  Select  <WINDOW>,  [USER  FILE 
NAME],  Place  cursor  at  the  top  of  the  first  open  column  -  right  of  “Verified” 
(Column  E). 

4.  CTRL  -  M;  to  start  analysis.  Analysis  is  completed. 

5.  A  Dialog  Box  appears  “POD  Curve  Table”.  Select  the  “maximum  crack  length  “  and 
“crack  increment”  to  be  used  in  plotting.  0.750  and  0.005  are  the  defaults. 

6.  Highlight  columns  $L  and  $M 

7.  From  the  <TOOL  BAR>  Select  [Insert  Chart} 

8.  A  Dialog  Box  Appears  -  Select  [xy  scatter];  NEXT 

9.  From  the  top  bar  Select  [Series]  NEXT . NEXT 

10.  A  dialog  box  appears.  -  Select  “Place  Chart”  [As  a  new  sheet], 

11.  Click  [Finish]-  POD  curve  appears. 

12.  From  the  <TOOL  BAR>  Select  Chart  -  Add  Data 

13.  Click  -  insert  new  data  -  a  dialog  box  appears  with  range. 

14.  Switch  to  sheet  1,  [Bottom  tool  bar]  select  $0:$P  -  [OK] 

15.  A  dialog  box  appears  {PASTE  SPECIAL] 

16.  Click,  x  “Categories  (X  values)  in  first  column). Then  OK 

17.  From  the  <TOOL  BAR>  Select  <CHART>  , CHART  TYPE> 

18.  A  dialog  box  appears.  Select,  “user  defined” 

19.  Click  previous  [USER  FILE  PLOT}  POD  —  OK.  (Must  have  loaded  the  format 
previously  -  To  load  the  applicable  mode-  See  Note  1 . 

20.  X’s  appear  on  the  top  and  bottom  of  the  chart  (Hit  /  Miss  Values).  Select  one  of  the 
X’s.;  Double  Click. 

21.  A  dialog  box  appears.  Select: 

22.  Line  should  be  NONE 

23.  Marker  should  be  CUSTOM 

24.  Select  Style  Symbol  (ie.  A) 

25.  Foreground  -  Select  BLACK 

26.  Background  -  Select  BLACK 

27.  Size  -  Select  8  pt. 

28.  [OK] 

29.  If  you  wish  to  have  an  insert  descriptive  box,  Select  from  another  chart  (be  sure  the 
outline  on  the  box  is  fully  cross  hatched.  Copy  and  paste  to  the  new  chart.  Text  can 
then  be  edited  to  annotate  the  new  chart. 
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~ .  POD 

Program  Instructions 
(Applicable  to  Excel  2000) 

30.  Enter  90%  POD  value  -  See  Note  2. 

31.  From  the  <TOOL  BAR>  Select  <FILE>  Click  Page  set-up,  Footer,  enter  description. 

32.  Data  analysis  is  complete. 

33.  Print  as  desired. 

Note  1: 

To  load  a  chart  template  (For  Example  -  Pen  1  A): 

a)  Load  a  previous  POD  plot  into  Excel. 

b)  From  Tool  Bar  >  Select  [CHART] 

c)  Select  [Chart  Type] 

d)  Select  [Custom  Type]  -  Top  Bar 

e)  Select  [User  Defined]—  Lower 

f)  Select  [Add] 

g)  Dialog  box  will  appear  -  INPUT  [chart  name  —  for  example  -  UT_Length 

h)  Select  -  [OK],  Dialog  Box  will  close 

i)  Use  this  user  format  in  subsequent  plotting. 


r\ 
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Note  2: 

To  Interpolate  for  the  90%  POD  Value 

a)  Got  to  Sheet  1  [  Data  ]  —  on  the  bottom  tool  bar. 

b)  Scroll  down  Column  M  to  the  values  at  89.XXX;  90.XXX 

c)  Select  these  four  cells  <Example> 


L 

M 

0.39 

89.86499 

0.395 

90.01687 

d)  Copy  these  four  cells. 

e)  From  the  <TOOL  BAR>  Select  <WINDOW>,  Select  interpol.xls 

f)  From  the  <TOOL  BAR>  Select  <EDIT>  “Paste  Special” 

g)  A  Dialog  box  appears.  Select  “Values” 

h)  The  interpolation  is  completed  and  the  90%  POD  value  appears  to  the  right  in  units  of 
both  inches  and  millimeters. 

i)  Copy  these  values  to  the  descriptive  data  box  on  the  <Chart>. 
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D.  COMMENTS 

The  Macro  is  relatively  simple,  but  was  designed  for  users  familiar  with  operation  of 
Excel  programs.  A  little  practice  is  required  to  become  comfortable  with  its  operation..  Application  to 
other  data  requires  modification  of  some  of  the  operating  parameters.  This  will  enable  analysis  of  data 
on  the  NDE  Capabilities  CD  is  practice  exercises. 

Electronic  copies  of  the  reference  paper  is  provided  on  the  CD  for  user  reference. 

REFERENCES: 

1 .  NDE  Capabilities  Data  Book  (3rd  Edition  -  1997),  available  from  NTIAC,  4 1 5  Crystal  Creek 
Drive,  Austin,  TX  78736-4725,  Telephone:  (512)  263-2106,(800)  NTIAC  39,  FAX:  (512)  263- 
3530;  CD  version  available. 

2.  Ward  D.  Rummel,  “Nondestructive  Evaluation  Engineering”,  A  Plenary  paper  prepared  for  and 
presented  as  an  “Indian  Society  for  Nondestructive  Testing,  Honorary  Fellow  Lecture  to  NDE  - 
99”,  17  December  1999,  Vadodara  (Baroda),  India. 
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